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^  This  report  investigates  the  validity  of  analytical  techniques  for  converting  continuous 
wave  (CW)  test  data  to  values  for  electromagnetic  pulse  (EMP)  to  measure  the  EMP  hard¬ 
ness  of  tactical  shelters.  MIL-STD-285  and  specifications  for  CW  testing  IEEE  299  specify 
measurements  in  the  near  field  using  dipole  and  loop  antennas.  Consequently,  the  test 
wavefronts  are  spherical,  rather  than  planar,  as  in  the  case  of  an  EMP.  Thus,  correction 
factors  for  EMP  needed  to  be  developed. 
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Two  methods  documented  in  the  literature  were  analyzed:  (1)  the  Schelkunoff  trans¬ 
mission  line  theory  approach  and  its  extension  to  slotted  shields,  and  (2)  the  plane  wave 
spectrum  approach  as  used  to  calculate  the  shielding  effectiveness  of  a  multiple-layer 
shield  to  the  fields  of  a  loop  antenna  positioned  in  the  coaxial  configuration.  It  was  found 
that  extension  of  the  transmission  line  theory  approach  to  slotted  shields  yielded  inaccurate 
results,  especially  for  loop  antennas.  However  the  plane  wave  spectrum  approach,  when 
combined  with  Huygen’s  principle,  yields  results  that  are  in  good  agreement  with  experi¬ 
mental  results. 

The  results  of  this  investigation  have  produced  recommendations  for  CW  test  techniques 
for  evaluating  EMP  hardness  of  tactical  shelters.  These  are: 

1.  Use  specifications  listed  in  IEEE  299,  with  modifications  as  given  in  the  Air  Force 
report  ESL-TR-80-01. 

2.  Make  maximum  use  of  the  IEEE  299  large-loop  test  or  seam  leak  detectors. 
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EXECUTIVE  SUMMARY 


Electromagnetic  fields  are  classified  as  magnetic,  plane,  or  electric  according  to  the  ratio 
of  their  magnetic  and  electric  field  components.  This  ratio  is  a  prime  factor  in  determining 
the  degree  to  which  a  field  will  be  reduced  by  an  electromagnetic  shield.  Within  the  electro* 
magnetic  pulse  (EMP)  spectrum  (100  Hz  to  100  MHz  for  a  high-altitude  nuclear  burst), 
it  is  easier  for  a  tactical  shelter  to  shield  against  a  plane  wave  field  than  a  magnetic  field; 
it  is  also  easier  to  shield  against  an  electric  field  than  a  plane  wave  field.  The  fields  caused 
by  a  high-altitude  nuclear  burst  are  plane  wave  fields.  The  standard  currently  used  to  verify 
the  EMP  hardness  of  tactical  shelters  (MIL-STD-285)  specifies  shielding  effectiveness 
measurements  with  magnetic  and  electric  fields  within  the  frequency  spectrum  of  an  EMP. 
Consequently,  shielding  effectiveness  measurements  made  on  a  tactical  shelter  according 
to  MIL-STD-285  will  only  give  upper  and  lower  bounds  on  the  shelter’s  EMP  shielding 
effectiveness.  The  spread  between  these  bounds  is  about  100  dB  in  the  EMP  spectrum. 
Furthermore,  the  upper  bound  provided  by  electric  field  measurements  is  often  greater 
than  the  dynamic  range  of  the  measuring  equipment. 

The  problems  associated  with  applying  MIL-STD-285  to  assessing  the  EMP  hardness 
of  shielded  enclosures  were  first  addressed  in  the  literature  by  Monroe  in  1973.*  Monroe 
applied  the  transmission  line  theory  of  shielding  and  developed  correction  factors  for 
MIL-STD-285  shielding  effectiveness  measurements  to  convert  them  to  values  meaningful 
for  EMP  plane  wave  fields.  The  only  other  literature  to  address  this  problem  (Villaseca  in 
1977**)  evaluated  the  validity  of  Monroe’s  work  from  an  analytical  standpoint.  Villaseca 
found  no  technical  errors  in  Monroe’s  work  within  the  stated  assumptions;  however,  he 
did  object  to  Monroe’s  application  of  transmission  line  theory  to  the  shielding  of  the  near 
fields  of  loop  and  monopole  antennas  (the  magnetic  and  electric  field  test  antennas,  respec¬ 
tively,  specified  in  MIL-STD-285).  Villaseca  derived  correction  factors  for  MIL-STD-285 
measurements  that  are  based  on  plane  wave  spectrum  theory  a  more  rigorous  analytical 
approach  than  transmission  line  theory.  These  correction  factors  agreed  with  Monroe's 
within  practical  limits.  Neither  Monroe  nor  Villaseca  experimentally  verified  their  cor¬ 
rection  factors. 

The  objective  of  this  research  was  to  apply  Monroe’s  and  Villaseca’s  analytical  tech¬ 
niques  to  problems  that  could  be  modeled  easily  in  CERL’s  electromagnetics  laboratory. 
Two  shielding  models  of  geometries  encountered  in  tactical  shelters  were  examined  : 
slotted  shields  and  multiple-layered  shields.  Computer  programs  were  developed  to 
predict  the  shielding  of  these  structures  for  the  fields  of  MIL-STD-285  test  antennas. 
A  computer  code  using  the  transmission  line  theory  outlined  in  Monroe’s  report  was  writ¬ 
ten  for  the  slotted  shield.  This  code  was  used  to  predict  the  shielding  that  was  measured 
on  two  slotted  shields  in  the  laboratory  with  loop  and  monopole  antennas.  Reasonable 
agreement  between  predicted  and  experimental  results  was  achieved  for  monopole  an¬ 
tennas;  however,  agreement  for  loop  antennas  was  not  acceptable.  Further  computer 
modeling  was  done  using  the  “Dipole-Moment  Approximation”  for  loop  antennas.  Im¬ 
proved  correlation  between  predicted  and  measured  results  was  obtained  for  loops  in 
both  the  coplanar  and  coaxial  test  orientations. 

Formulations  based  on  the  plane  wave  spectrum  approach  were  used  in  a  computer 
program  to  calculate  the  shielding  effectiveness  of  multiple-layered  shields  for  the  fields 
of  a  loop  antenna  in  the  coaxial  test  configuration.  The  shielding  predicted  by  the  program 


*R.  L.  Monroe,  EMP  Shielding  Effectiveness  end  MIL-STD-285,  HDL-TR-1636/AD771997 
(Harry  Diamond  Laboratories.  U.S.  Army  MaterielCommand,  July  1973). 

**E.  Villaseca,  C.  Davis,  W.  Blackwood,  and  W.  Getson,  An  Investigation  of  the  Validity  of  Applying 
MIL-STD-285  to  EMP  Shielding  Effectiveness,  ADA051889,  prepared  by  Harris  Corp.  Electronics 
System  Division  (Defense  Nuclear  Agency,  15  April  1977). 


was  compared  with  measurements  taken  in  the  laboratory  on  dual-layer,  copper-foil 
shields  mounted  in  an  aperture  in  a  steel-construction  shielded  room.  When  the  effects 
of  this  aperture  were  accounted  for  with  a  separate  computer  program,  agreement  with 
experimental  results  was  excellent. 

Direct  experimental  evaluation  of  the  applicability  of  MIL-STD-285  to  the  assessment 
of  the  EMP  hardness  of  tactical  shelters  could  not  be  made  in  this  study.  However,  at  this 
time,  the  best  techniques  for  testing  a  tactical  shelter’s  shielding  effectiveness  are  applica¬ 
tion  of  the  procedures  outlined  in  MIL-STD-285  and  in  IEEE  299,  along  with  the  guide¬ 
lines  presented  in  this  report.  Particular  attention  must  be  paid  to  electromagnetic  leakage 
at  discontinuities  in  the  shield  panels,  such  as  seams,  door  seals  and  cable  entries.  Also, 
test  frequencies  outside  the  EMP  spectrum  may  need  to  be  specified,  depending  on  the 
mission  of  the  shelter. 
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EVALUATION  OF  APPLICABILITY  OF 
STANDARD  CW  EMI/RFI  SHIELDING 
EFFECTIVENESS  TEST  TECHNIQUES  TO 
ASSESSMENT  OF  EMP  HARDNESS  OF 
TACTICAL  SHELTERS 


1  INTRODUCTION 


Background 

The  accepted  technique  for  determining  the  electro¬ 
magnetic  pulse  (EMP)  hardness  of  a  shielded  enclosure 
or  tactical  shelter  is  to  subject  it  to  a  threat-level  pulse 
environment.  This  is  generally  not  a  practical  solution, 
however,  because  EMP  simulation  equipment  is  expen¬ 
sive.  MIL-STD-2851  and  IEEE  299 2  list  the  radio 
frequency,  continuous  wave  (CW)  test  techniques 
agreed  on  by  shielding  manufacturers  and  by  personnel 
who  accept  shielded  enclosures.  The  measurements 
specified  by  these  standards  are  made  in  the  near  field 
using  dipole  and  loop  antennas.  Consequently,  the  test 
wavefronts  impinging  on  the  shielded  enclosure  are 
spherical,  rather  than  planar,  as  in  the  case  of  an  EMP. 
Specifically,  the  wave  impedance  of  the  test  wavefronts 
will  be  above  and  below,  respectively,  that  of  the  EMP 
planewave  for  the  dipole  and  loop  fields.  Thus,  the 
shielding  effectiveness  measured  by  the  dipole  and  loop 
tests  will  be  above  and  below  that  of  an  incident  EMP. 

Monroe3  has  proposed  extension  of  the  Shelkunoff 
transmission  line  theory  and  Villaseca4  has  proposed  use 
of  the  plane  wave  spectrum  approach  in  an  attempt  to 
develop  correction  factors  for  converting  ’  ‘IL-STD-285 
measurements  to  values  meaningful  for  EMP  planar 
radiation.  However,  knowledge  of  these  efforts  is  not 


'Military  Standard  Attenuation  Measurements  for  Enclo¬ 
sures,  Electromagnetic  Shielding  for  Electronic  Test  Purposes, 
Method  of,  MIL-STD-285  (Department  of  Defense,  25  June 
1956). 

2 Proposed  IEEE  Recommended  Practice  for  Measurement 
of  Shielding  Effectiveness  of  High-Performance  Shielding 
Enclosures,  IEEE  299  (Institute  of  Electrical  and  Electronics 
Engineers  (IEEE1,  June  1969). 

3R.  L.  Monroe,  EMP  Shielding  Effectiveness  and  MIL- 
STD-285,  HDL-TR-1636/AD771997  (Harry  Diamond  Labora¬ 
tories,  U.S.  Army  Materiel  Command.  July  1973). 

*E.  Villaseca,  C.  Davis,  W.  Blackwood,  and  W.  Getson, 
An  Investigation  of  the  Validity  of  Applying  MIL-STD-285  to 
EMP  Shielding  Effectiveness,  ADA051889,  prepared  by  Harris 
Corp.  Electronics  System  Division  (Defense  Nuclear  Agency, 
15  April  1977). 


widespread,  and  the  results  have  not  been  applied  to 
the  procurement  of  shielded  facilities.  Furthermore, 
little  experimental  work  has  been  done  to  verify  the 
correction  factors. 

The  military  currently  has  no  CW  test  techniques  to 
directly  measure  the  EMP  hardness  of  tactical  shelters; 
furthermore,  there  has  been  no  experimental  verifica¬ 
tion  of  previous  efforts  to  convert  erusting  accepted 
standard  CW  measurements  to  values  meaningful  for 
EMP.  Therefore,  the  U.S.  Army  Construction  Engineer¬ 
ing  Research  Laboratory  (CERL)  was  requested  to 
examine  the  validity  of  the  analytical  techniques  used 
to  convert  CW  test  data  to  values  for  EMP. 

Objective 

The  objectives  of  this  study  were  to;  (l)examine  the 
validity  of  the  EMPcorrection  factors  for  MIL-STD-285 
measurements  developed  by  Monroe  and  Villaseca,  and 
(2)  recommend  CW  test  procedures  for  verifying  EMP 
hardness  of  tactical  shelters. 

Approach 

A  literature  search  was  conducted  to  provide  informa¬ 
tion  on  previous  work  done  to  provide  EMP  correction 
factors  for  MIL-STD-285  measurements  (Chapter  2). 
Computer  programs  were  used  to  calculate  the  shield¬ 
ing  effectiveness  predicted  by  analytical  models  for 
slotted  shields  and  multiple-layer  shields.  Laboratory 
experiments  were  then  carried  out  to  verify  the  results  of 
these  computer  programs  (Chapters  3  and  4).  The  effects 
of  cavity  resonances  on  tactical  shelters  were  investi¬ 
gated  (Chapter  5).  CW  test  procedures  for  verifying 
EMP  hardness  of  tactical  shelters  were  recommended, 
based  on  the  results  of  the  research  (Chapter  6). 

Mode  of  Technology  Transfer 

It  is  anticipated  that  this  study  will  impact  on  the 
application  of  the  radio  frequency  interference  (RF1) 
shielding  effectiveness  measurement  techniques  speci¬ 
fied  in  MIL-STD-285.  IEEE  299  and  NSA  65-6  for 
the  EMP  hardness  testing  of  GLCM  control  electronics 
modules. 


2  SUMMARY  OF  CORRECTION  FACTORS 
DEVELOPED  FOR  EMP 


Cost  constraints  and  the  nonavailability  of  portable 
EMP  sim,  lators  often  make  full-scale,  threat-level  EMP 
hardness  testing  of  shielded  facilities  and  tactical  shel- 
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ters  impossible.  Usually,  a  threat-level  pulse  test  would 
be  considered  only  for  initial  acceptance  or  design 
evaluation  testing.  CW  shielding  effectiveness  test  tech¬ 
niques  offer  economical  alternatives  to  pulse  tests  for 
EMP  hardness  testing  of  production  shelter  units  and 
for  periodic  maintenance  testing. 


B  =  20  log 
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MIL- STD-285  specifies  the  test  procedures  most 
often  used  to  measure  the  electromagnetic  shielding 
effectiveness  of  shielded  enclosures.  IEEE  299  outlines 
similar  procedures  with  additional  tests  for  seam  leakage 
(large  loop  tests).  Within  the  generally  accepted  fre¬ 
quency  spectrum  of  an  EMP  caused  by  a  high-altitude 
nuclear  burst  (100  Hz  to  100  MHz),  the  procedures 
specitied  in  M1L-STD-285  and  'EEE  299  call  for  use 
of  loop  and  monopole  antennas  as  radiating  sources  for 
the  measurements.  These  give  shielding  effectiveness 
measurements  for  low-impedance  (magnetic)  and  high- 
impedance  (electric)  fields,  respectively.  The  wave 
impedance  of  a  high-altitude  EMP  impinging  on  a  tacti¬ 
cal  shelter  will  be  377  J2,  which  lies  between  the  wave 
impedances  of  the  loop  and  monopole  sources  at  the 
specified  measurement  distances,  and  within  the  EMP 
spectrum.  Thus,  the  question  is  what  information  the 
standard  test  procedures  of  MIL-STD-285  and  IEEE 
299  can  give  about  the  EMP  shielding  effectiveness  of 
tactical  shelters. 

The  problem  of  applying  MIL-STD-285  test  proce¬ 
dures  to  evaluating  the  EMP  shielding  effectiveness  of 
shielded  enclosures  was  first  considered  by  Monroe 
in  EMP  Shielding  Effectiveness  and  MIL-STD-2855. 
Monroe  employs  Schelkunoff  s  transmission  line  theory 
approach  to  shielding  to  address  the  fact  that  the  mis¬ 
match  between  the  intrinsic  impedance  of  the  shielding 
material  and  the  wave  impedance  of  the  impinging 
radiation  is  critical  to  a  structure's  shielding  effective¬ 
ness.  This  fact  is  apparent  in  the  standard  transmission 
line  theory  equation  for  shielding  effectiveness: 

SE  =  R  +  A  +  B  (dB)  [Eq  1  j 


where: 


R  =  20  log 


ik  +  n* 

4|k| 


.  initial  reflection  loss 


a  =  (nyof )%,  reciprocal  of  shield  skin  depth 


f  =  frequency 
t  =  shield  thickness 
M  =  shield  permeability 
o  =  shield  conductivity 

i=y/-T 

The  term  corresponding  to  the  initial  reflection  loss,  R, 
depends  directly  on  the  ratio  of  the  impinging  wave 
impedance  to  the  intrinsic  impedance  of  the  shield. 

Monroe  observes  that  MIL-STD-285  measurements 
made  with  loops  and  monopoles  give  lower  and  upper 
bounds,  respectively,  on  a  shielding  enclosure's  EMP 
shielding  effectiveness  at  a  given  frequency.  This  is 
because  the  impedance  mismatch  between  a  loop  source, 
an  EMP  (plane  wave),  a  monopole  (or  dipole)  source, 
and  the  intrinsic  impedance  of  a  shield  is  ordered  as: 


^loop 

S' 

Zemp 

S' 

^monopole 

ZshWd 

ZjhHd 

V* 

Zshidd 

The  spread  between  the  theoretical  bounds  given  by  loop 
and  monopole  measurements  can  be  about  100  dB  in 
the  frequency  region  betow  1  MHz.  Thus,  a  structure's 
actual  EMP  shielding  effectiveness  cannot  be  accurately 
estimated  directly  from  MIL-STD-285  measurements 
alone.  (Furthermore,  neither  MIL-STD-285  or  IEEE 
299  specify  both  loop  and  monopole  tests  at  any  given 
frequency.) 


A  =  8.686  art,  atte.iuation  loss  in  penetrating  shield 
once 


5R.  L.  Monroe,  EMP  Shielding  Effectiveness  and  MIL- 
STD-285,  HDL-TR-t 6  36/AD77 1 997  (Harry  Diamond  Labora¬ 
tories,  U  S.  Army  Materiel  Command,  July  1973). 


Monroe  has  evaluated  the  transmission  line  expres¬ 
sions  for  the  difference  oetween  the  shielding  of  a 
given  conductoi  to  an  EMP  and  to  a  loop  source 
(or,  equivalently  between  an  EMP  and  a  monopole 
[dipole]  source);  he  then  gives  a  correction  factor 
5(0  which  can  be  added  to  a  loop  measurement  or  sub- 
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traded  from  a  monopole  measurement  to  estimate 
the  shelter's  EMP  shielding  effectiveness: 


where:  Zemp  —  377  Q 


Thus,  6(0  can  be  evaluated  at  the  test  frequency  for  the 
appropriate  source  and  added  10  the  MIL-STD-285 
test  result  in  order  to  estimate  the  structure’s  actual 
EMP  shielding  effectiveness. 


reactive,  nonpropagating  energy  stored  in  the  near  field 
of  the  antenna.  The  transmission  line  approach  used  by 
Monroe  does  not  allow  consideration  of  the  reactive 
energy  which,  in  tne  near  field,  is  relatively  large  com¬ 
pared  to  the  real,  propagating,  energy.  Viliaseca  states 
that  the  reactive  coupling  of  energy  is  of  fundamental 
importance  in  the  problems  associated  with  shielding 
effectiveness  measurements  made  in  the  near  field  with 
low-level  CW  sources. 

In  the  plane  wave  spectrum  approach,  an  electro¬ 
magnetic  field  (E  or  H)  is  expanded  in  a  set  of  plane 
waves  as: 

(Eq  4) 


The  validity  of  Monroe’s  work  is  addressed  by 
Viliaseca.6  This  study  acknowledges  that  Monroe's 
work  was  the  first  step  in  applying  standard  EMI 
shielding  effectiveness  measurement  techniques  to  the 
EMP  problem.  It  concludes  that,  within  the  stated 
assumptions,  there  are  no  technical  errors  in  the  solu¬ 
tion  presented.  However,  Viliaseca  disagrees  with  the 
underlying  assumption  on  which  the  report  is  based. 
Monroe’s  application  of  the  transmission  line  theory 
of  shielding  inherently  assumes  that  a  plane  wave  is 
incident  upon  the  shielding  material.  However,  Viliaseca 
points  out  that  this  assumption  is  not  true  for  the 
actual  situation  in  the  MIL-STD-285  test  setup,  since 
the  shield  is  in  the  extreme  near  field  of  t!  j  test  an¬ 
tennas.  When  the  source  antenna  is  placed  close  to  the 
shield,  the  incident  waves  are  of  a  complex  spherical 
nature,  rather  than  planar. 


E(x.  y,  z)  =  JT  E(KX,  Ky)e-j(K*x+Kyy+K*2>dKx  dKy 
where: 

E  (Kx,  Ky)  =  the  plane  wave  spectrum  or  density  of 
the  electric  field  per  square  wave  number 

K  =  Kx  X  +  Ky  9  +  K.z  2  is  the  propagation 
vector. 


From  Maxwell’s  equations,  it  follows  that 

2k 

®-'"x 

which  implies  that 


[Eq  5) 

[Eq  6] 


Viliaseca  derives  correction  factors  for  MIL-STD-285 
measurements  via  an  independent  technique,  free  of 
any  simplifying  assumptions;  this  technique  uses  the 
concept  of  an  ’’Angular  Spectrum  of  Plane  Waves.” 
Introduced  by  Booker  and  Clemmou7  in  1957,  the 
concept  allows  the  representation  of  an  antenna’s 
radiation  pattern  by  a  linear  combination  of  plane 
waves.  The  steady-state  radiation  pattern  of  the  antenna 
reduces  to  an  angular  spectrum  of  these  plane  waves 
taken  over  all  angles.  The  portions  of  the  pattern  func¬ 
tion  that  correspond  to  complex  angles  determine  the 

*E.  Viliaseca,  C.  Davis,  W.  Blackwood,  and  W.  Getson, 
An  Investigation  of  the  Validity  of  Applying  MIL-STD-285  to 
EMP  Shielding  Effectiveness,  ADA0S1889,  prepared  by  Harris 
Corp.  Electronics  System  Division  (Defense  Nuclear  Agency, 
15  April  1977). 

tH.  G.  Booker  and  P.  L  Clemmou,  Proceedings  of  IEEE, 
Vol  97,  Part  HI  (1957),  pp  11-17. 


As  shown  by  Eq  6,  the  (Kx,  Ky)  plane  can  be  divided 
into  two  regions  separated  by  a  circle  centered  at  the 
origin  of  radius  0.  The  inner  part  of  the  circle,  corre¬ 
sponding  to  real  values  of  Kz,  defines  the  propagating 
portion  of  the  plane  wave  spectrum  and  is  called  the 
“visible”  region.  The  outer  part  of  the  circle,  called  the 
“invisible”  region,  corresponds  to  imaginary  values  of 
Kz  and  defines  the  non-propagating,  evanescent  waves 
which  are  exponentially  attenuated  for  z>0. 

In  short,  Villaseca's  approach  derives  transmission 
coefficients  for  the  visible  and  invisible  portions  of 
transverse  magnetic  (TM)  and  transverse  electric  (TE) 
plane  wave  power  spectrum  expansions  for  dipole  and 
loop  test  antennas,  respectively.  Then,  the  shielding 
effectiveness  of  a  plane  shield  is  determined  by  calcu¬ 
lating  the  ratio  of  the  power  transmitted  through  the 
shield  to  the  power  incident  on  the  shield  in  both  the 
visible  and  invisible  regions  of  the  spectrum.  The  same 
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calculations  are  also  performed  for  a  uniform  plane 
wave  incident  on  the  shield  in  order  to  produce  correc¬ 
tion  factors  applicable  to  MIL-STD-285  measurements. 
Villaseca  found  that  if  the  calculations  were  confined 
to  the  visible  region,  correction  factors  agreed  precisely 
with  those  of  Monroe.  Including  the  contribution  of  the 
power  in  the  invisible  region,  produced  a  loop  correc¬ 
tion  factor  about  10  dB  greater  and  a  dipole  correction 
factor  about  10  dB  less  than  Monroe's. 

The  plane  wave  spectrum  technique  provided  correc¬ 
tion  factors  that  did  not  agree  precisely  with  those 
derived  by  the  transmission  line  technique;  however, 
Villaseca  concluded  that  Monroe’s  report  represented  a 
good  engineering  approximation  for  the  case  of  a  uni¬ 
form  conductivity  shield  with  no  apertures  or  seams.  He 
cautions,  however,  against  the  application  of  Monroe’s 
approach  to  shields  with  penetrations  or  seams. 

In  summary,  two  basic  approaches  have  been  used  to 
obtain  correction  factors  for  converting  MIL-STD-285 
loop  and  monopole  (dipole)  shielding  effectiveness 
measurements  to  values  meaningful  for  an  incident 
EMP  (plane  waves).  The  Schelkunoff  transmission 
line  approach  used  by  Monroe  represents  a  first-order 
approach  to  the  problem  and  is  considered  to  be  a 
reasonable  approximation  only  for  a  uniform  planar 
shield.  As  shown  in  Chapter  3,  the  extension  of 
Monroe’s  approach  to  slotted  shields  does  not  always 
yield  accurate  results. 

The  plane  wave  spectrum  concept  used  by  Villaseca 
is  a  more  rigorous  and  versatile  analytical  approach  for 
calculating  shielding  effectiveness.  Since  it  does  not  rely 
on  any  simplifying  assumption,  it  holds  great  promise 
for  application  to  general  classes  of  shielding  problems. 
Chapter  4  discusses  the  application  of  plane  wave  spec¬ 
trum  techniques  to  calculating  the  shielding  effective¬ 
ness  of  double-layered  shields  when  illuminated  by  a 
loop  antenna  in  the  coaxial  orientation.  Appendix  A  pre¬ 
sents  an  application  of  plane  wave  spectrum  techniques 
to  calculating  the  shielding  effectiveness  of  a  multiple- 
layer  shield  when  illuminated  by  arbitrary  fields. 


3  SHIELDING  EFFECTIVENESS  OF 
SLOTTED  SHIELDS 


slot  is  one  whose  length  is  much  larger  than  its  width. 
This  discussion  is  limited  to  slots  that  are  electrically 
small;  i.e.,  one  whose  length  is  much  smaller  than  the 
wavelengths  of  the  impinging  radiation. 

Monroe8  extends  the  Schelkunoff  transmission  line 
theory  of  shielding  to  the  calculation  of  the  shielding 
effectiveness  of  a  perfectly  conducting  shield  contain¬ 
ing  an  electrically  small,  narrow  slot.  Based  on  the 
assumption  that  the  illuminating  radiation  will  be 
approximately  uniform,  Monroe  states  (without  further 
justification)  that  the  shielding  effectiveness  due  to 
reflection  from  the  slot  can  be  expressed  as: 

lk+llJ 

SE  =  20,Og'W  [Eq?1 

where:  k  = 

^sk>t 

and  Z,|ot  =  impedance  of  the  slot. 

This  equation  has  exactly  the  same  form  as  the  reflec¬ 
tion  term  in  the  expression  for  the  shielding  effectiveness 
of  a  planar  shield  with  no  discontinuities.  For  a  slotted 
shield,  the  intrinsic  impedance  of  the  shielding  material 
has  been  replaced  by  the  slot  impedance. 

The  slot  impedance  is  highly  dependent  on  the  polar¬ 
ization  of  the  incident  field.  Maximum  coupling  through 
the  slot  is  obtained  when  the  incident  electric  field  vec¬ 
tor  is  aligned  in  a  direction  perpendicular  to  the  large 
dimension  of  the  slot.  For  this  polarization,  the  slot 
impedance  is  related  to  the  driving  point  impedance, 
Z^j,  of  the  complementary  dipole  and  is  given  by 
Monroe  as: 


•7  _  V72  _  T}2  Red  "  j^cd 

islo'  =  4Zcd=4 


(Eq81 


where: 

Tj=  377a 

Red  =  real  part  of  Z^ 

Xed  *  imaginary  part  of  Zjd 


The  narrow  slot  is  a  common  analytical  model  for 
discontinuities  in  practical  shielded  enclosures  such  as 
welds  and  door  seams.  In  the  present  context,  a  narrow 


*R.  L.  Monroe,  EMP  ShMdl/it  Effectirtnttt  titd  MIL- 
STD-285,  HDL-TR-1 6  36/AD771 997  (Harry  Diamond  Labora¬ 
tories.  u.S.  Army  Materiel  Command.  July  1973). 
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Approximate  expressions  for  R^j  and  Xcd  are  cited  by 
Monroe  as: 


Z0  sinh  (2y) 

^  Z  cosh2  (7)  -  cos2  (0L) 

[Eq9| 

Z0  sin  (20L) 

^cd  Z  cosh2  (y)  -  cos2  (0L) 

[Eq  10] 

where: 

kHz  to  300  MHz.  To  compare  the  predictions  of  the 
transmission  line  theory  to  the  experimental  data. 
Eqs  7  through  10  were  programmed  on  a  CDC  7600 
computer  via  the  FORTRAN  program  SLOT  (see 
Appendix  B).  Shielding  effectiveness  values  for  loop  and 
monopole  sources  located  1  ft  away  from  the  shield  were 
calculated  for  both  slot  sizes.  Experimental  data  were 
taken  with  the  loop  antennas  oriented  in  the  coplanar 
position  with  the  plane  of  the  antennas  perpendicular 
to  the  slot’s  long  dimension.  Data  for  the  monopoles 
were  taken  with  the  antennas  oriented  perpendicular  to 
the  slot.  These  antenna  orientations  provided  an  inci¬ 
dent  field  with  the  electric  field  vector  oriented  in  a 
direction  perpendicular  to  the  slot’s  long  dimension. 


An  examination  of  the  behavior  of  the  slot  impedance 
with  frequency  indicates  that  within  the  EMP  spectrum, 
the  shielding  effectiveness  of  a  slotted  shield  is  approxi¬ 
mately  constant  with  frequency  for  the  fields  produced 
by  a  loop  source.  The  shielding  of  the  fields  produced 
by  a  dipole  or  monopole  by  a  slotted  shield  shows  about 
a  40-db/decade  decrease  with  increasing  frequency. 

CW  testing  was  conducted  on  slotted  shield  panels 
in  the  laboratory  to  examine  the  accuracy  of  the  trans¬ 
mission  line  theory  of  shielding  when  applied  to  slots. 
Two  1/4-in.  thick  aluminum  panels,  one  with  a  1-meter 
by  1/8-in.  slot  and  the  other  with  a  1 /2-meter  by 
1/ 16-in.  slot,  were  mounted  over  a  2-  x  4-ft  test  aperture 
in  a  steel  shielded  room.  The  panels  were  tested  in  the 
MIL- STD-285  configuration  with  I2-in.  loopand41-in. 
monopole  antennas  over  the  frequency  band  from  100 


Figures  1  and  2  present  the  experimental  results, 
along  with  the  transmission  line  theory  predictions,  for 
the  shielding  effectiveness  of  the  slotted  panels  to  the 
monopole  fields  for  both  slot  sizes.  As  shown  in  the 
figures,  the  transmission  line  theory  predictions  are 
within  an  average  of  6  dB  of  the  measured  values,  for 
the  larger  slot,  while  the  predictions  for  the  smaller 
slot  show  about  a  14  dB  underestimate  of  shielding 
effectiveness. 

Figures  3  and  4  compare  the  transmission  line  theory 
predictions  of  shielding  effectiveness  for  loop  antennas 
in  the  coplanar  orientation  and  show  experimental 
measurements.  For  this  source  and  orientation,  the 
transmission  line  theory  predictions  underestimate  the 
shielding  effectiveness  of  the  slotted  panel  by  about  14 
and  18  dB  for  the  larger  and  smaller  slots,  respectively. 
Thus,  for  loops  in  the  coplanar  orientation  (a  common 
method  for  checking  seam  leakage),  the  slot  impedance 
shielding  model  does  not  appear  to  give  accurate  results. 

C.  M.  Butler9  suggests  an  alternate  method  for 
calculating  the  shielding  effectiveness  of  a  perfectly 
conducting  shield  with  an  electrically  small,  narrow 
slot.  The  technique,  known  as  the  “Dipole-Moment 
Approximation.”  allows  the  expression  of  the  electro¬ 
magnetic  fields  that  penetrate  an  electrically  small 
aperture  in  terms  of  the  fields  of  equivalent  electric 


9C.  M.  Butler,  “Dipole  Moment  Approximation  and 
Polarizabilities,”  Section  2.1. 3. 2  in  EMP  Interection  Principitt. 
Technique!,  end  Reference  Dele,  EMP  Interaction  2-1,  AFWL- 
TR-80-402  (Air  Force  Weapons  Lab,  Air  Force  Systems 
Command,  December  1980). 
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SHIELDING  EFFECTIVENESS  («)  SHCLDING  EFFECTIVENESS  «•) 
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Figure  1.  Shielding  effectiveness  of  slotted  shield  to  monopole  fields.  Slot  size  =  1  m  X  1/8  in. 


FREQUENCY 

Figure  2.  Shielding  effectiveness  of  slotted  shield  to  monopole  fields.  Slot  size  *  1/2  m  X  1/16  in. 
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SHIELDING  EFFECTIVENESS  (46) 


60 


O 


SHIELDING  EFFECTIVENESS  OF  SLOTTED  SHIELD 
TO  LOOP  FIELDS,  COPLANAR  ORIENTATION 

SLOT  SIZE  *  I  METER  X  1/6  INCH 
MEASURED  VALUES  O 

PREDICTED  VALUES  A 


o 


o 


o 


o 


10 


A 


lOOKHz 


I  MHz 


10  MHz 


FREQUENCY 


OOMA 


Figure  3.  Shielding  effectiveness  of  slotted  shield  to  loop  fields,  coplanar  orientation.  Slot  size  llmX  1/8  in. 
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Figure  4.  Shielding  effectiveness  of  slotted  shield  to  loop  fields,  coplanar  orientation.  Slot  size  *  1/2  m  X  1/16  in 


Figure  5.  Orientation  of  slot  in  cartesian  coordinate 
system. 


netic  polarizability  tensors,  at  and  atm  ,  by: 

Pa  =  2  tte  •  Ejc>  m,  =  -2am  •  H*  [Eq  12j 
where: 

Ejj.,  Hk  are  the  short  circuit  fields  in  the  aperture 
(equal  to  two  times  incident  fields) 

S;  and  are  two-dimensional  symmetric  dyads. 

The  components  of  the  polarizability  tensors  for  a 
narrow  slot  oriented  in  a  cartesian  coordinate  system 
as  shown  in  Figure  5  are  given  on  page  437  of  Butler.10 
Calculations  for  the  slot  dimensions  used  in  this  study 
yield 

*x  =  2-132  X  l<r2»  1.979  X  1(T*  = 

at,  zz  ~  ttjn,  yy 

for  the  1 -meter  by  1/8-in.  slot  and 

ocm «  =  2.665  X  1CTJ  »  2.474  X  10"1  = 

ae,  zz  =  am.  yy 


and  magnetic  dipoles  in  the  aperture.  These  dipoles 
can  be  expressed  in  terms  of  the  fields  incident  on  the 
aperture.  The  remainder  of  this  discussion  will  be 
confined  to  the  shielding  of  the  fields  of  loop  antennas 
that  are  positioned  in  both  the  coplanar  and  coaxial 
orientations.  Consequently,  only  the  magnetic  field 
propagating  from  the  aperture  will  be  considered. 

The  magnetic  field,  H*1,  radiated  by  an  electric  dipole, 
p,,  and  i  magnetic  dipole,  in,,  located  at  the  origin  of  a 
cartesian  coordinate  system  in  free  space,  is  given  by: 

Hd(r,s)=-s  p,(s)X  VG(r,s)-VX  [Eq  11] 
(in,  (s)  X  V  G  (r,  s)] 

where. 

G  (r ,  s)  =  e_i  T  r/4trr 
X  =  wavelength 

The  aperture  equivalent  electric  and  magnetic  dipoles, 
pa  and  ma.  are  given  in  terms  of  the  electric  and  mag- 


for  the  1/2 -meter  by  1/1 6-in.  slot.  Thus,  the  contribu¬ 
tion  to  the  magnetic  field  Hd  by  the  electric  dipole  is 
negligible,  and  Eq  1 1  reduces  to 

H°(r,  s)=-VX  [ma (s) X  V  G (r,  s)]  (Eq  13] 

where  only  the  om  xx  component  of  crm  is  non-zero. 
Thus,  ma  is  given  by 

in,  =  -4  Om,  xx  Hjnc,  x  *  [Eq  14] 

where: 

Hj„c,  x  the  x-component  of  the  magnetic  field  inci¬ 
dent  on  the  aperture.  Letting 

*  -4  am,  xx  Hmc.  x  (Eq  15] 


10C.  M.  Butler,  “Dipole  Moment  Approximation  and 
Polaris*  bill  tie*,"  Section  2.1. 3.2  in  BMR  Interaction  Mnctptet, 
Technique *,  and  Reference  Data,  EMP  Interaction  2-1,  AFWL- 
TR-80-402  (Air  Force  Weapons  Lab,  Air  Force  Systems 
Command,  December  1980). 


and  substituting  into  Eq  7  yields 


loops  via  the  standard  Field  equations  for  an  elementary 
loop.  Shielding  effectiveness  was  then  calculated  for 
the  coplanar  orientation  as: 


e-jkr  e-jkr\ 

+  -rrr  x 


H**  (f)  (m*9y*  4ffr  T  “‘X3x2  4rrr  ) 

(  32 

+  \  m*3x3y  4rrr  /  ^ 


/  J2_  e-jkr 

\  m*3xdz  4irr 


■* 


[Eq  16] 


Referring  to  the  coordinate  system  shown  in  Figure  5, 
the  components  of  If*  of  interest  are  Hj  for  the  coaxial 
orientation  and  H£  for  the  coplanar  orientation. 


The  derivatives  indicated  in  Eq  16  were  calculated 
using  a  quasi-static  approximation  for  the  exponential 
term.  The  result  for  the  x-component  of  Hd  (coplanar 
loops)  is 


H--' 


m«  e"*kr 

3(yJ  +  zJ) 

4ir 

(x2  +  y2  +  z2)*72 

[Eq  17] 


The  result  for  the  z-component  (coaxial  loops)  is 


* 


mx  e~lkr  3xz 
4rr  (x2  +  y2  +  z2)v2 


[Eq  18] 


To  calculate  these  magnetic  field  components,  Eqs  17 
and  18  were  integrated  over  the  slits  after  making  the 
substitutions 


x  =  x  -  x' 

y  =  y  -  y' 

z  =  z  -  z’  [Eq  19] 

where  the  primed  coordinate  system  has  its  origin  in 
the  aperture  and  the  unprimed  system  is  centered  at 
the  observation  point  (the  center  of  the  receiving 
loop).  These  integrals  were  carried  out  over  the  primed 
variables  via  standard  numerical  techniques  in  two 
FORTRAN  programs  SLITCP  (SLIT  CoPlanar)  and 
SLJTCA  (SLIT  CoAxial).  (Listings  of  SLITCP  and 
SLITCA  appear  in  Appendix  B.)  Reference  field 
strengths  were  calculated  at  the  center  of  the  receiving 


jjRcf 

SE  =  20  log  -j-jj- 


and  for  the  coaxial  orientation  as: 

jjRef 

SE  =  20  log 


[Eq  20] 


[Eq  21] 


Figures  6  and  7  present  the  results  of  the  program 
SLITCP  and  the  experimental  measurements  taken  for 
the  coplanar  loops.  Also  included  are  the  predictions 
of  the  transmission  line  theory  approach  from  Figures 
3  and  4.  The  agreement  between  the  Dipole-Moment 
Approximation  technique  and  the  experimental  mea¬ 
surements  approaches  1  dB  for  the  smaller  slot  and 
represents  about  a  17-dB  improvement  in  accuracy 
over  the  transmission  line  theory  approach. 


Figures  8  and  9  given  the  measurements  taken  with 
loop  antbnnas  in  the  coaxial  orientation  and  the 
predictions  of  the  computer  program  SLITCA.  For 
this  orientation,  the  agreement  between  experimental 
results  and  the  Dipole-Moment  Approximation  predic¬ 
tion  again  improves  slightly  as  the  slot  size  decreases. 
For  the  smaller  slot,  the  average  of  the  experimental 
results  is  46  dB,  while  SLITCA  predicts  42.6  dB  of 
shielding  over  the  frequency  band. 

In  summary,  the  results  of  these  experimental  mea¬ 
surements  showed  that  extension  of  the  Shelkunoff 
transmission  line  theory  of  shielding  to  slotted  shields 
yielded  inaccurate  results,  especially  for  coplanar  test 
loops.  A  more  promising  technique  is  the  Dipole- 
Moment  Approximation,  whichyieldedaccuratepredic- 
tions  for  electrically  small  slots.  Butler11  has  extended 
the  Dipole-Moment  Approximation  and  related  tech¬ 
niques  to  deal  with  more  complicated,  practical  open¬ 
ings,  such  as  hatch  apertures. 


"C.  M.  Butter,  “Dipole  Moment  Approximation  and 
Polarizabilities,’’  Section  2. 1.3. 2  in  EMF  Interaction  Mndpiet, 
Techniques  and  Reference  Data,  EMP  Interaction  2-1,  AFWL- 
TR-S 0-402  (Air  Force  Weapons  Lab,  Air  Force  Systems 
Command,  December  1980). 
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Figure  7.  Shielding  effectiveness  of  slotted  shield  to  loop  fields,  coplanar  orientation.  Slot  size  *  1/2  m  X  1/16  in. 
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Figure  9.  Shielding  effectiveness  of  slotted  shield  to  loop  fields,  coaxial  orientation.  Slot  size  *  1/2  m  X  1/16  in 


4  SHIELDING  EFFECTIVENESS  OF 
MULTIPLE-LAYER  SHIELDS 


The  construction  of  the  Ground  Launch  Cruise 
Missile  shelter  outer  casing  is  a  composite  of  plane 
sheets  of  laminated  conductor  and  structural  materials. 
The  conductor  sheets  provide  electromagnetic  shielding 
while  at  least  one  filler  layer  is  used  to  increase  rigidity. 
There  is  a  lack  of  research  literature  about  the  develop¬ 
ment  of  the  capability  to  calculate  the  shielding  effec¬ 
tiveness  of  plane-laminated  materials.  This  chapter 
describes  a  technique  developed  during  this  investiga¬ 
tion  for  calculating  the  shielding  effectiveness  of  mul¬ 
tiple-layer  shields  for  loop  antennas  in  the  coaxial 
configuration.  This  technique  is  based  on  the  plane-wave 
spectrum  approach.  Also  included  is  a  discussion  of  the 
computer  programming  of  these  relations,  as  well  as 
experimental  assessment  of  their  validity  and  accuracy. 

Tyras12  gives  the  transfer  function  for  plane  wave 
propagation  through  multiple  planar  layers  of  arbitrary 
material.  Figure  10  shows  the  geometry  and  notation 


12  George  Tyras,  Radiation  and  Propagation  of  Electro¬ 
magnetic  Waves  (Academic  Press,  1969),  pp  40-43. 
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Figure  10.  Geometry  and  notation  for  multiple-layer 
planar  shield. 


used.  The  transfer  function  of  the  multiple  plane  layers 
is  given  by: 


—  _  2  Zmtl  m 

*  *  n  .  *■»  «* 


Vi 


S,  +  Z,  p=2  Sp  cos  ♦p  -  i  Zp  sin  *p 

c  -7  /  Sp  -  ‘  tan  \ 
Vl  %,-iS.  tan*J 


•iZptan^p 

-  i  Sp  tan  *p  , 
*P  =  VdP 


[Eq  22] 

[Eq  23] 
[Eq  24] 


where: 


p  =  material  layer:  p=l,  2...,  m+1 
p=l  is  the  incident  half-space 
p=m+!  is  the  transmission  half-space 

m+1  =  total  number  of  layers  (including  half-spaces) 

m-1  =  total  number  of  shielding  and  structural  layers 

dp  =  depth  of  a  shielding  or  structural  layer 

kpz  =  complex  propagation  constant  in  the  z  direction 
for  layer  p 

Zp  *  intrinsic  impedance  of  a  layer 

Sp  =  input  impedance  at  p*  interface 

Eq  23  is  a  recursive  relation  for  calculating  Sp.  The  cal¬ 
culation  begins  by  setting  the  input  impedance  at  the 
last  interface  equal  to  the  intrinsic  impedance  of  the 
transmission  half-space: 


V=Zn 


[Eq  25] 


Computational  effort  is  reduced  by  combining  Eqs 
22  and  23  to  obtain: 


T  = : 


m 

n 


1 


1  +  ?,  p=2  cos  ♦p  (1  +  i  £p  tan  *p) 
in  which 


[Eq  26] 


r. 


p  z„ 


[Eq  27] 


The  values  for  fp  are  obtained  by  dividing  Eq  23 
through  by  Zp.,  and  rearranging  to  obtain: 


fp-i 


-  hzl ,  _?E_  *p  -  j tin  V 


Zp-i  Zp_ j  1  -  i  fp  tan  ♦p 


[Eq  28] 


20 


in  which  the  recursion  is  now  started  by  setting 


a  =  (5.8  X  107)or,  conductivity  of  the  shield 


ft-,  -  7*^-1  fEq  29J 


at  =  relative  conductivity  of  the  shield  with 
respect  to  copper 


The  transfer  function  given  in  Eq  26  applies  to  plane 
waves  incident  upon  a  multiple-layer  medium.  It  is 
extendable  to  the  case  of  the  nonplanar  radiation  of 
a  loop  antenna  if  the  fields  of  the  loop  antenna  are 
expressed  as  an  angular  spectrum  of  plane  waves. 


Moser13  provides  the  derivation  of  the  shielding 
effectiveness  of  a  single-layer,  planar  shield  for  loop 
antennas  in  the  coaxial  configuration.  This  work  was 
later  pursued  by  Bannister.14  The  equation  for  shield¬ 
ing  effectiveness,  SE,  is: 


4u 

SE  »  -20  log  ~(a3  +  z3)*3. 

rc\2T 

/  ~T  Jt(Xa)exp{-T<>z-t(r  -  t0) }dX 

Jo  To 


[Eq  30) 


a  =  radius  of  the  loop  antenna 

z  =  ri  +  r2>  the  center-to-center  separation  of 
the  transmitting  and  receiving  loop  antenna 

Jt(Xa)  =  the  Bessel  function  of  order  one  and 
argument  Xa. 

The  Bessel  function  term  corresponds  to  a  plane  wave 
spectrum  description  of  the  fields  due  to  the  loop  source. 

Eq  3  may  be  expressed  in  terms  of  the  transfer  func¬ 
tion  for  the  shield,  T(X): 


SE  *  -20  log 


(a3  +  z3)*3 
a 


where: 


C  =  [(t/t„  +  Hr)1  -  (t/t„  -  Hr)1  e_2tT]_l 

r  =  (X3  +  r3)1/3 

To  =  (X+y1o)l/i 


T(X)  J ,  (Xa)exp  {-r0(z-t) }  dX  [Eq  3 1  ] 


where: 


T(X) 


4^Ct  -t. 


[Eq  32] 


7o  =  i  2w/Xlir,  the  free  space  propagation 
constant 

=(iwii0/iro)1/2,  the  propagation  constant  in 
the  shield 

t,  =  |or|/{2  =(w/uoAir0/2)l/3  =  1/5 
5  -(2/v/p0pro)in,  the  skin  depth  in  the  shield 
p ,  =  relative  permeability  of  the  shield 
t  *  shield  thickness 

H0  =  4wX  10~7,  permeability  of  free  space 


l3I.  R.  Moser,  “Low-Frequency  Shielding  of  e  Circular 
Loop  Electromagnetic  Field  Source,”  IEEE  Trenmctiom  on 
Electromagnetic  Compatibility,  Vol  EMC-9,  No.  1  (March  1969). 

,4P.  R.  Bannister,  “Further  Notes  for  Predicting  Shielding 
Effectiveness  for  the  Plane  Shield  Case,”  IEEE  Trvmctiom  on 
Electromagnetic  Compatibility,  Vol  EMC-11,  No.  2  (May  1969). 


After  some  algebraic  operations,  it  may  be  shown  that 
the  multiple-layer  transfer  function  given  by  Eq  26 
reduces  to  Eq  32  for  the  case  m+1  *  three  total  layers. 
This  is  the  case  for  a  single  plane  shield,  with  air  as 
the  half-space  layers  on  each  side.  Eq  31  is  also  the 
generalized  version  of  the  expression  for  the  shielding 
effectiveness  of  a  single  layer  shield;  this  equation  is 
useful  for  computing  the  shielding  effectiveness  of  the 
multiple  layer  case  once  the  transfer  function  T(X)  is 
calculated  according  to  Eq  26. 

Before  Eq  26  can  be  used  in  Eq  31.  one  must  convert 
the  notation  for  Tyres'  propagation  constant  to  that 
used  by  Moser  and  Bannister.  Tyres  defines  the  square 
of  the  propagation  constant  as: 

ka  *  <jj2fie  +  ioup  [Eq  33 1 

while  Bannister  uses  the  definition 

y3  =  -o o3fie  +  jow»i  [Eq  34J 

These  two  notations  can  be  related  by  letting 
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v  =  -k5 

[Eq  35] 

j  =-i 

[Eq  36] 

The  change  indicated  in  Eq  36  is  dictated  by  the  fact 
that  the  time  convention  used  in  this  work  is  ei“t,  as 

opposed  to  e'kJt  as  used  by  Tyras. 

Now  the  multiple-layer  transfer  function  given  by 
Eq  26  must  be  expressed  in  terms  of  the  angular  fre¬ 
quency  integration  constant,  X,  used  in  Eq  31.  First, 

'Pp  in  Eq  24  is  written  as  a  function  of  rp . 

'Pp  ~  kpzdp  i  Tpdp  -  j  rpdp 

[Eq  37] 

Next,  rp  is  writ  -/i  as  a  function  of  X: 

rp  -  VXJ  +  y2 

[Eq  38] 

Putting  Eqs  37  and  38  into  Eq  26  and  converting  to 
hyperbolic  functions  yields: 


To  assess  the  validity  and  accuracy  of  MULSH, 
shielding  effectiveness  measurements  were  taken  in 
the  laboratory  on  panels  composed  of  two  1-mil-thick 
copper  sheets  separated  by  either  a  1/8-in.  or  1/2-in. 
thickness  of  masonite  (nonconductive  material  for 
structural  support).  These  panels  were  mounted  over 
the  2-  x  4- ft  aperture  in  a  shielded  room.  LP-105  loop 
antennas  (12-in.  diameter)  were  used  for  measurements 
greater  than  100  kHz;  22  1/2-in.  diameter  loops  made 
at  CERL  were  used  for  measurements  at  and  below 
100  kHz. 

Initial  comparisons  between  experimental  shielding 
effectiveness  measurements  and  MULSH  predictions 
indicated  that  MULSH  consistently  underestimated 
shielding  effectiveness  in  the  10-kHz  to  1-MHz  range 
by  about  10  dB.  However,  this  was  not  surprising,  since 
MULSH  models  an  infinitely  large  (height  and  width) 
double-copper- foil  panel,  whereas  the  experimental 
setup  consisted  essentially  of  a  2-  x  4-ft,  double-foil 
aperture  in  an  otherwise  opaque  screen. 


[Eq  39] 


T(\)  = 


m 

n 


1  +  f,  p=2  cosh(rpdp)[l  +  fptanh  (rpdp)] 
{•„  +  tanh  (r„d„) 


=  l  [Eq  4 1  ] 


where 


To  evaluate  MULSH  more  completely,  another 
program  was  developed  to  account  for  the  geometry  of 
the  aperture  in  the  experimental  setup.  This  program, 
called  APRAD  (APerture  RADiation),  calculates  the 
radiation  of  fields  of  a  loop  antenna  (coaxial  configura¬ 
tion)  through  an  aperture  containing  a  single  thin  con¬ 
ducting  sheet  in  an  otherwise  opaque  screen.  APRAD  is 
listed  in  Appendix  B.  Although  APRAD  was  developed 
for  this  special  case,  the  theoretical  techniques  used  in 
its  derivation  can  be  extended  to  the  general  problem 
of  calculating  radiation  through  apertures. 


_  Mrp  rp-i 

Zp-i  Mr(p-I)  rP 


[Eq  42] 


has  been  substituted  into  Eq  40.  This  is  the  form  of  the 
multiple-layer  transfer  function  used  in  the  integral 
expression  (Eq  31). 


The  above  equations  were  programmed  on  a  CDC 
7600  computer  via  the  FORTRAN  program  MULSH 
(MULtiple  SHieldsXsee  Appendix  B).  Using  dX  =  0.10, 
the  integral  expression  in  Eq  31  converged  to  a  relative 
error  of  1C4  within  1000  steps.  (This  control  parameter 
may  have  to  be  changed  if  problem  dimensions  differ 
appreciably  from  those  studied  in  this  work.)  MULSH 
allows  the  user  to  specify  a  shield  design  having  one  to 
five  total  shielding  and/or  structural  layers  of  any  thick¬ 
ness  and  material.  (It  should  be  noted  that  MULSH 
assumes  an  infinitely  large  shield.) 


The  formulas  used  to  develop  APRAD  are  based  on 
the  application  of  Huygen’s  principle.15  The  fields  on  the 
transmission  side  of  the  foil  aperture  can  be  expressed 
in  terms  of  the  tangential  electric  field  impressed  on 
the  incidence  side  of  the  foil  by  the  transmitting  loop 
antenna.  Since  the  receiving  loop  responds  to  the 
magnetic  field  normal  to  the  plane  of  the  loop,  only 
the  magnetic  field  is  calculated  on  the  transmission  side 
of  the  screen.  This  is  given  by: 

H(r)  =  O^o)"'  ( V  X  V  X  ?(?))  [Eq  43] 


where: 


i5R.  F.  Harrington,  Time-Harmonic  Electromagnetic  Fleldt 
(McGraw-Hill  Book  Company,  Inc.,  1961),  p  110. 
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and 


M  (r')  =  T(X)  (2E,*,,  X  n);  magnetic  current  in  foil 
to  =  2irf;  f  =  frequency  in  Hz 
k  =  2ir/X 

e0  =  8.854  X  10*u  F/m 

r'  =  (x'x  +  y'y):  vector  from  origin  in  plane 
of  aperture  to  a  point  in  the  plane  of  the 
aperture 

r  =  (xx  +  yy  +  z):  vector  from  origin  in 

plane  of  aperture  to  a  point  in  the  plane  of 
the  receiving  antenna  located  R,*.  meters 
from  aperture 

T(X)  =  transfer  function  of  sheet  as  calculated 
in  Eq  39 

X  =  v/u2  +  v2:  angular  frequency  variable 
u  =  (2»r/X)  sind  cos  <t> 
v  *  (2tr/X)  sind  sin  4> 

d,  =  angles  of  a  spherical  coordinate  system 

with  origin  at  the  center  of  the  transmitting 
loop  and  Z  axis  coaxial  with  the  transmitting 
and  receiving  loops. 

As  outlined  by  the  above  equations,  the  approach 
taken  in  APRAD  is  as  follows.  First,  the  tangential 
electric  field  is  calculated  on  the  incidence  side  of  the 
aperture  from  the  standard  equations  for  the  fields 
produced  by  a  small  (compared  to  wavelength)  loop. 
Then  the  electric  field  is  passed  through  the  foil  using 
the  transmission  coefficient  of  Eq  39  for  a  single  layer; 
the  coefficient  is  evaluated  at  values  of  X  corresponding 
to  angles  that  the  aperture  intercepts.  In  the  case  of  the 
2-  x  4-ft  aperture,  the  electric  field  is  evaluated  on  a 
5x9  sample  point  grid  centered  on  the  aperture.  The 
transfer  function  is  the  same  one  used  in  MULSH; 
however,  in  APRAD,  it  is  evaluated  only  for  angular 
frequencies  corresponding  to  angles  intercepted  by 
the  aperture.  In  MULSH,  an  infinite  sheet  is  assumed. 
Thus,  the  transfer  function  in  MULSH  appears  as  an 
integral  over  X,  with  an  infinite  limit  of  integration. 

After  the  electric  field  is  passed  through  the  foil, 
a  magnetic  current  flowing  in  the  foil  is  expressed  as 


a  cross  product  of  the  reduced  electric  field  and  the 
normal  vector  to  the  foil  surface.  Then  an  electric  vec¬ 
tor  potential,  F  (dual  of  a  magnetic  vector  potential), 
is  expressed  as  an  integral  (sum)  over  the  sample  points 
of  the  magnetic  current  in  the  foil.  From  the  electric 
vector  potential,  the  magnetic  field  in  the  plane  of  the 
receiving  loop  antenna  is  calculated  via  a  double  curl 
operation. 

In  performing  the  derivatives  for  the  double  curl, 
a  quasi-static  approximation  was  made  by  assuming 
that  the  exponential  term  within  the  integral  for  F  was 
constant.  This  was  justified,  since  |r- r*|  varies  by  less 
than  a  factor  of  2  over  the  aperture,  and  k  has  a  maxi¬ 
mum  value  of  .021  in  the  frequency  band  of  interest.  The 
normal  component  of  the  magnetic  field  is  calculated 
within  the  receiving  loop  at  7  points  and  integrated  to 
give  the  received  signal  strength.  To  calculate  shielding 
effectiveness,  a  reference  response  is  calculated  for  the 
loop  antenna  (separated  by  an  appropriate  distance 
of  free  space)  directly  from  the  field  equations  for  a 
small  loop.  The  normal  component  of  the  magnetic 
field  is  again  calculated  at  7  points  in  the  receiving  loop 
and  integrated. 

To  correct  the  MULSH  predictions  for  the  experi¬ 
mental  setup,  APRAD  and  MULSH  are  both  run  for  a 
single  1-mil  copper  sheet.  The  difference  between  the 
shielding  effectiveness  predicted  by  the  two  programs 
is  termed  the  “Blockage  Factor”  and  accounts  for 
the  opaque  part  of  the  shielded  room  wall.  Computa¬ 
tional  results  have  shown  that  the  Blockage  Factor 
for  this  experimental  setup  is  about  10  dB,  thus 
accounting  for  MULSH’s  10-dB  underestimate  of 
shielding  effectiveness. 

Figure  1 1  shows  the  predicted  values  of  shielding 
effectiveness  produced  by  MULSH  and  APRAD  for  a 
single,  1-mil-thick,  copper  sheet.  Also  shown  are  experi¬ 
mental  measurements  for  the  single  copper  sheet  and 
the  difference  between  the  predictions  of  MULSH  and 
APRAD  (the  Blockage  Factor).  Note  that  the  APRAD 
predictions  show  reasonable  correlation  with  the  experi¬ 
mental  measurements  except  at  10  MHz,  where  the 
experimental  accuracy  is  highly  suspect. 

Figure  1 2  shows  the  MULSH  predictions  as  corrected 
by  addition  of  the  Blockage  Factor  for  the  two  double¬ 
copper  foil  shields.  Notice  the  increase  in  shielding 
effectiveness  as  the  separation  between  the  sheets  is 
increased  from  1/8-in.  to  1/2  in.  This  can  be  explained 
by  examining  the  behavior  of  S,(Eq  23).  the  input  impe¬ 
dance  at  the  first  air-copper  interface.  As  the  separation 
between  the  copper  sheets  increases.  St  decreases 
(numerical  results  from  MULSH).  Consequently,  the 
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Figure  11.  Shielding  effectiveness  of  single  copper  sheet  (1-mil-thick)  to  loop  fields,  coaxial  orientation. 
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Figure  12.  Shielding  effectiveness  of  double  copper  shield  to  loop  fields,  coaxial  orientation. 
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mismatch  between  the  impedance  of  the  multiple 
layer  shield  and  the  impedance  of  the  incident  wave  is 
increased.  This  increases  the  reflection  from  the  inter¬ 
face,  thus  increasing  shielding  effectiveness. 

In  summary,  a  plane  wave  spectrum  approach  was 
used  to  calculate  the  shielding  effectiveness  of  a  mul¬ 
tiple  layer  shield  to  the  fields  of  a  loop  antenna  posi¬ 
tioned  in  the  coaxial  test  configuration.  To  model  the 
aperture  of  the  experimental  setup  more  accurately, 
Huygen’s  principle  was  used  to  express  the  fields  on 
the  transmission  side  of  the  aperture  (containing  a 
single,  thin-layer  shield)  in  terms  of  the  fields  incident 
on  the  aperture.  When  these  two  approaches  are  com¬ 
bined  (APRAD  and  MULSH,  and  the  Blockage  Factor), 
the  results  are  in  good  agreement  with  the  experimental 
results.  Appendix  A  presents  an  extension  of  the  plane 
wave  spectrum  approach  to  calculating  the  shielding  of 
a  multiple-layer  medium  to  arbitrary  incident  radiation. 


5  CAVITY  RESONANCE  EFFECTS 


A  screen’s  EMP  shielding  effectiveness  can  often  be 
estimated  from  its  CW  performance.  However,  for  a 
cavity-type  structure  (e.g.,  an  enclosed  shelter),  be¬ 
havior  near  its  interior  resonances  cannot  be  directly 
extrapolated  to  estimate  the  strength  of  EMP  penetra¬ 
tion  into  the  shelter.  The  effect  of  resonance  on  the 
CW  field  penetration  into  the  interior  of  cavity-type 
configurations  has  been  investigated  and  documented 
by  Mittra14  and  by  Yung.17 

Figure  13  shows  that  the  field  at  the  center  of  an 
infinitely  long  cylinder  can  become  quite  large  near  its 
interior  resonance.  Figure  14,  which  show'  the  induced 
current  in  a  wire  located  inside  a  cylindrical  cavity  with 
an  aperture,  also  exhibits  a  similar  phenomenon.  How¬ 
ever,  to  compute  the  transient  response  of  these  struc¬ 
tures,  one  must  first  multiply  the  frequency  response 
of  the  structure  and  the  spectrum  of  the  incident  EMP 


14R.  Mittra,  “Cavity  Excitation  via  Apertures,”  EMP  Inter¬ 
action:  Principles,  Techniques  and  Reference  Data.  Ch.  2. 3.2.1, 
K.  S.  H.  Lee,  ed.  (Air  Force  Systems  Command,  December 
1980),  pp  522-523. 

17E.  K.  Yung,  S.  W.  Lee,  and  R.  Mittra,  “Penetration  of  an 
EM  Wave  into  a  Cylindrical  Cavity  and  the  Current  Induced  on 
a  Wire  inside,”  AEU.  Band  33.  Heft  4  (1979),  pp  149-156. 


field.  Butler1*  provides  an  example  ol  such  a  compu¬ 
tation.  Figure  15  shows  the  frequency  response  of  a 
screen  with  an  aperture;  Figure  16  exhibits  the  transient 
response  of  the  same  response  for  an  incident  EMP 
whose  waveform  is  shown  in  the  inset  of  the  same 
figure.  As  shown,  the  primary  effect  of  the  resonance 
is  to  introduce  ringing  in  the  transient  response;  the 
frequency  of  the  ringing  is  determined  by  the  resonant 
frequency  of  the  cavity,  while  the  decay  envelope  is 
determined  by  the  width  of  the  resonance.  The  energy 
deposited  by  the  incident  pulse  within  the  resonant 
bandwidth  of  the  cavity  primarily  determines  the  peak 
of  the  ringing  waveform.  Consequently,  the  height  of 
this  peak  is  not  enhanced  nearly  as  dramatically  as  the 
resonance  peak  of  the  frequency  response  curve;  this  is 
because  the  incident  EMP  pulse  seldom  contains  much 
energy  in  the  resonance  band  of  the  cavity. 

The  worst  case  energy  coupling  into  a  tactical  shelter 
at  a  cavity  resonance  frequency  can  be  roughly  calcu¬ 
lated  as  follows.  For  typical  shelter  dimensions,  the 
lowest  cavity  resonance  will  be  around  100  MHz.  Figure 
17  illustrates  the  normalized  power  spectrum  for  an 
EMP;  the  energy  content  at  100  MHz  is  60  dB  below 
the  peak  value  of  .075  V/m-liz  at  700  kHz.  Thus,  the 
power  density  incident  on  the  shelter  in  a  1-Hz  band 
width  is  given  by: 

„  1  (.075  X  10~3)3 

d  2  377 

=  7.46  X  10“ 12  w/m3  -  Hz  [Eq  44) 

If  the  EMP  is  incident  on  a  shelter  with  a  wall  of,  say, 
3  x  7  m.  then,  with  no  shielding,  the  power  coupled 
into  the  shelter  would  be. 

P  =  (21  m3X7.46  X  10'13  w/m3  -  Hz) 

=  1.57  X  10'10w/Hz  [Eq  45) 

for  a  1-Hz  band  width.  Even  for  a  1-kHz  band  width, 
this  is  still  only  .157  jrW  of  power  in  the  resonant 
band,  without  considering  the  shielding.  Since  a  high- 
performance  shelter  will  probably  have  at  least  100  dB 
of  shielding  to  a  plane  wave  at  100  MHz,  cavity  reso¬ 
nances  of  tactical  shelters  are  not  expected  to  be  stimu¬ 
lated  by  EMPs. 


'*C.  M.  Butler,  Y.  Rahmat-Samii.  and  R.  Mittra,  “Electro¬ 
magnetic  Penetration  Through  Apertures  in  Conducting  Sur¬ 
faces,”  IEEE  Trans,  on  Antennas  and  Propagation,  Voi  AP-26, 
No,  1  (January  1978),  pp  82-93. 
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Figure  13.  Variation  of  interior  field  along  0  =  0°,  with  ka  as  parameter  for  field  at  the  center  of  the  cavity  versus  ka 
for  0O=  10°  and  30°. 
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Figure  14.  Induced  current  at  the  center  of  the  wire 
inside  the  cavity  as  a  function  of  cavity 
length  2h.The  geometrical  parameters  are: 
a  =0.5  X,  c  =  0.3  X,  d  =  0.015  X, 
a*,  =  0.001  X,  h*  =  0.2  X,  e/a  =0.1. 
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Figure  15.  Frequency  domain  behavior  of  |E^  (f.  r_)| 
sampled  at  a  point  2  m  behind  the  aperture. 
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6  SHELTER  TEST  PROCEDURES 
RECOMMENDATIONS 


Discussion 

As  stated  in  Chapters  1  and  2,  when  it  is  not  possible 
or  practical  to  use  an  EMP  simulator  to  evaluate  shielded 
structures,  rooms,  zones,  etc.,  it  is  customary  to  sub¬ 
stitute  test  methods  described  in  MIL-STD-285,19 
IEEE  299,30  or  NSA  65-6.31 .  These  and  similar  pro¬ 
cedures  require  testing  at  discrete  CW  frequencies 
using  small  antennas  such  as  loops,  dipoles,  monopoles 
and  horns.  Problems  or  shortcomings  of  this  type  of 
testing  include  the  following: 

1.  Using  loops  or  dipoles  for  CW  testing  involves 
testing  with  either  low  (loop)  or  high  (dipole  or  mono¬ 
pole)  wave  impedance  fields  through  most  of  the  fre¬ 
quency  range  of  interest.  Monroe”  and  Villaseca33  have 
shown  that  shielding  effectiveness  values  resulting  from 
these  measurements  can  be  related  to  shielding  effec¬ 
tiveness  for  plane  waves  where  the  wave  impedance  is 
3  77  ohms.  These  conversions,  however,  have  not  gained 
wide  acceptance. 

2.  The  EMP  transient  has  spectral  energy  distributed 
over  a  relatively  broad  frequency  range.  For  complete 
testing,  a  continuous  frequency  scan  should  be  taken 
of  the  entire  EMP  frequency  range.  Experience  has 
shown,  however,  that  the  shielding  effectiveness  of 
shielded  structures  such  as  rooms  and  shelters  tends  to 
follow  a  relatively  smooth  curve  at  frequencies  below 
20  MHz.  Therefore,  the  curve  up  to  20  MHz  can  be 

,9MUitary  Standard  Attenuation  Measurements  for  Enclo¬ 
sures,  Electromagnetic  Shielding  for  Electronic  Test  Purposes. 
Method  of,  MIL-STD-285  (Department  of  Defense,  25  June 
1956). 

30 Proposed  IEEE  Recommended  Practice  for  Measurement 
of  Shielding  Effectiveness  of  High-Performance  Shielding 
Enclosures,  IEEE  299  (Institute  of  Electrical  and  Electronics 
Engineers  [IEEEj,  June  1969). 

3. P.  R.  Trybus,  national  Security  Agency  Specification  for 
R.F.  Shielded  Enclosures  for  Communications  Equipment: 
General  Specification,  NSA  65-6  (National  Security  Agency, 
30  October  1964). 

33R.  L.  Monroe,  EMP  Shielding  Effectiveness  and  MIL- 
STD-285,  HDL-TR-1636/AD771997  (Harry  Diamond  Labora¬ 
tories,  U.S.  Army  Materiel  Command,  July  1973). 

3JE.  Villaseca,  C.  Davit,  W.  Blackwood,  and  W.  Getton, 
An  Investigation  of  the  Validity  of  Applying  MIL-STD-285  to 
EMP  Shielding  Effectiveness,  ADA051889,  prepared  by  Harris 
Corp.  Electronics  System  Division  (Defense  Nuclear  Agency, 
15  April  1977). 


approximated  by  obtaining  a  few  measurements  at 
discrete  frequencies.  However,  at  frequencies  above 
20  MHz,  problems  in  taking  the  CW  shielding  measure¬ 
ments  begin  to  occur.  These  problems  are  related  to 
reflections  and  standing  waves  within  the  shielded 
volume  being  tested,  which  cause  large  variations  and 
lack  of  repeatability  of  test  data.  The  reflections  are 
further  complicated  by  the  effects  of  test  personnel 
movements  and  the  changes  in  measured  test  data 
caused  by  small  changes  in  transmitting  antenna  orien¬ 
tation  relative  to  a  defect.  Thus,  trying  to  obtain  a 
good,  average  CW  shielding  effectiveness  value  above 
20  MHz  is  time-consuming  and  requires  skilled  and 
persistent  test  operators.  At  best,  the  data  taken  above 
20  MHz  will  have  doubtful  repeatability  within  5  dB 
and  even  worse  repeatability  as  frequency  is  increased. 

3.  Another  problem  in  CW  testing  is  the  inability  to 
account  for  the  effect  of  numerous  cavity  resonances. 
These  resonances  are  defined  by  the  mathematical 
expression: 

ie,“' 

where: 

a,  b,  and  c  are  interior  dimensions  of  the  enclosure 
in  meters. 

m,  n,  and  p  are  integers,  only  one  of  which  may  be 
zero  at  a  time. 

For  typical  room-sized  enclosures,  the  lowest  resonance 
is  near  100  MHz.  but  resonances  are  lower  for  larger 
enclosures.  The  problems  caused  by  the  resonances 
are:  (a)  apparent  shielding  effectiveness  at  a  specific 
resonance  may  be  much  lower  than  actual  shielding 
effectiveness  at  frequencies  away  from  resonances, 
and  (b)  it  is  difficult  to  use  CW  test  techniques  for  all 
resonances  without  a  sweeping/tracking  transmitter/ 
receiver  system.  In  an  EMP  simulator,  the  transient  pulse 
has  spectral  components  which  may  excite  all  enclosure 
resonances  simultaneously;  thus,  the  combined  reso¬ 
nance  effects  may  be  measured  for  each  simulator  pulse. 
However,  simulator  tests  cannot  generally  provide  the 
absolute  wont-case  resonance  effects  either,  because 
the  amount  of  excitation  of  wave  guide  modes  within 
the  enclosure  is  uncertain. 

Summary  of  the  State  of  the  Art  in  CW  Testing 

Consideration  of  the  current  state  of  the  art  in  CW 
shielding  effectiveness  testing  should  account  for  (1) 
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standard  test  specifications  and  recent  work  to  upgrade 
them,  (2)  test  methods  regarded  within  the  industry 
as  “acceptable,”  and  (3)  recent  work  to  develop  new 
approaches.  These  are  discussed  below. 

Standard  Test  Specifications 

Currently,  there  are  two  standard  test  specifications 
(MIL-STD-285  and  NS  A  65-6)  and  one  proposed 
standard  (IEEE  299).  No  modifications  to  these  stan¬ 
dards  have  been  published  since  their  initial  printings. 
The  former  U.S.  Army  Electronics  Command  did 
some  preliminary  work  in  the  mid-1970s  to  revise 
MIL-STD-285,  but  it  was  not  completed.  Other 
efforts  to  revise  the  IEEE  299  standard  are  as  yet 
unpublished.14  A  preliminary  draft  of  these  revisions 
has  noted  that  no  significant  changes  to  the  testing 
methods  are  recommended.  There  have  been  no  known 
attempts  to  revise  NSA  65-6. 

Recently,  the  Mission  Research  Corporation  prepared 
a  “Specification  for  EMP  Attenuating  Prototype  Air 
Force  Tactical  Shelter”15  for  the  Air  Force  Weapons 
Laboratory  (AFWL).  This  specification,  which  includes 
a  section  on  EMP  testing,  states  that  three  separate 
tests  should  be  done  on  the  shelter:  (1)  extended 
MIL-STD-285  testing,  (2)  EMP  simulation  tests,  and 
(3)  EMP  hardness  maintenance  checks.  The  extended 
MIL-STD-285  tests  differ  from  MIL-STD-285  only  in 
that  more  frequencies  are  specified  in  H-field,  E-field, 
and  plane  wave  testing,  and  more  test  point  locations 
are  specified.  The  EMP  simulation  tests  require  a  EMP 
simulator.  The  EMP  hardness  maintenance  checks  are 
done  with  a  radio  frequency  (RF)  seam  sniffer. 

Specifications  for  CW  testing  of  SAFEGUARD  facil¬ 
ities  have  been  published  by  the  Army. 16  These  methods 
use  the  same  principles  set  forth  in  MIL-STD-285  with 
some  slight  modifications. 

Test  Methods  Used  by  Industry 

Shielded  room  manufacturers  have  used  some  sig¬ 
nificantly  abbreviated  test  methods  for  EMI  shielding 

^Proposed  Revision  of  IEEE  Recommended  Practice  for 
Measurement  of  Shielding  Effectiveness  of  High  Performance 
Shielding  Enclosures  (Draft)  (IEEE). 

15  A  G.  Finei,  H.  M.  Fowles,  P.  R.  Trybus,  Ground  Based  C’ 
Facilities  Hardening  and  Validation,  GBC*-3-MRC-064 
(Minion  Research  Corporation,  November  1980). 

2*H.  E.  Atkins,  R.  E.  Evans,  B.  1.  Gay,  H.  L.  Holt,  and 
A.  R.  Wright,  Safeguard  Tactical  Ground  Facilities  EMP/RFI 
Facilities  Acceptance  Construction /Installation  Test  Plan- 
Grand  Forks,  prepared  by  the  Boeing  Company  (U.S.  Army 
Engineer  Division,  Huntsville,  6  January  1972). 


evaluations.  They  feel  that  these  abbreviated  procedures 
indicate  EMP  shielding  adequately.  The  procedures 
use  a  seam  leak  detector  to  locate  leaks  and  establish 
relative  signal  leakage  magnitudes.  The  seam  leak 
detection  method  uses  CW  current  injection  at  about 
100  kHz;  the  current  injection  points  are  located 
on  diagonally  opposite  comers.  Generally,  one  set  of 
measurements  is  taken  with  injection  on  one  set  of 
comers,  but  repeated  using  another  set.  Wont  readings 
are  taken  as  representative. 

Typically,  shelter  manufacturers  closely  follow  the 
procedures  of  MIL-STD-285  to  measure  EMI /R FI 
shielding  effectiveness.  However,  they  use  the  mini¬ 
mum  number  of  antenna  test  positions  allowed  by 
MIL-STD-285  around  the  shelter. 

CW  testing  has  been  used  to  verify  EMP  hardness  of 
large  facilities,  such  as  Minuteman  Missile  structures, 
the  SAFEGUARD  ABM  structures,  the  Systems  Tech¬ 
nology  Test  Facility  (STTF)  on  Meek  Island,  and 
numerous  secure  communications  facilities. 

Development  of  New  CW  Test  Methods  , 

Two  programs  are  considering  CW  testing  for  EMP 
hardness  evaluation.  AFWL  is  developing  specialized 
antennas  called  plane  wave  launchers.  AFWL  has  already 
developed  a  test  concept  known  as  the  PARTES  con¬ 
cept,  in  which  a  planar  array  of  small  antenna  elements 
(loop  or  monopole)  generates  a  propagating  wave  with 
a  nearly  planar  wavefront.  These  antenna  arrays  can 
be  used  with  CW  excitation  and  will  allow  CW  testing. 
This  can  provide  actual  plane  wave  shielding  effective¬ 
ness  data  without  the  uncertainty  caused  by  a  non- 
planar  (spherical)  wavefront,  as  generated  by  loop  and 
dipole  antennas.  However,  the  problem  of  accounting 
for  resonances  is  still  not  solved  by  the  plane  wave 
test  approach. 

In  another  program,  the  Defense  Nuclear  Agency 
(DNA)  is  investipting  the  use  of  a  sophisticated  trans¬ 
mitter  and  receiver  system  which  has  a  very  high 
measurement  range  and  the  capability  to  obtain  swept 
frequency  measurements  between  10  kHz  and  100 
MHz.  This  system  can  measure  shielding  effectiveness 
magnitudes  according  to  MIL-STD-285  specifications 
and  also  measure  the  phase  function  over  the  entire 
frequency  range  in  discrete  frequency  steps  smaller 
than  1  Hz.  Thus,  it  is  possible  to  obtain  magnitude  and 
phase  curves  for  CW  shielding  effectiveness  over  the 
frequency  ranp  containing  more  than  99.9  percent 
of  the  energy  of  a  typical  high-altitude  EMP  (HEMP) 
transient.  Fourier  analysis  can  be  used  to  predict  the 


EMP  transient  response  from  these  curves,  obviously 
excluding  the  effects  of  resonances  and  shielding  effec¬ 
tiveness  above  100  MHz. 

In  both  of  these  programs,  the  systems  are  in  the 
developmental  phases  and  not  ready  for  recommenda¬ 
tion  for  tactical  shelter  maintenance  testing.  Reviews 
of  the  CW  test  methods  are  available  in  the  literature.3’ 

Summary 

The  following  discussion  summarizes  the  factors 
that  led  to  test  procedure  recommendations  discussed 
on  pp  30  through  31 . 

The  equipment  needed  for  the  new  systems  and 
concepts  being  studied  by  AFWL  and  DNA  is  not 
readily  available  for  the  widespread  use  required  for 
periodic  maintenance  testing  of  tactical  shelters.  Fur¬ 
thermore,  these  systems  will  probably  be  too  expensive 
for  the  intended  applications  and  will  require  a  greater 
amount  of  operator  skill  and/or  data  analysis  than  is 
desirable  for  maintenance  testing. 

In  most  applications  for  shelters  of  International 
Standards  Organization  (ISO)  size  or  smaller,  the  effects 
of  cavity  resonances  should  not  be  significant,  because 
more  than  99.9  percent  of  the  HEMP  energy  is  below 
the  lowest  resonant  frequency.  (See  Chapter  5.) 

Experience  has  shown  that  at  frequencies  below  20 
MHz,  the  shielding  effectiveness  versus  frequency  curve 
is  relatively  smooth  and  without  discontinuities.  Thus, 
it  can  be  approximated  by  measurements  at  a  few 
(three  or  four)  appropriately  spaced  frequencies. 

The  shielding  effectiveness  measured  for  the  CW 
nonplane  wave  can  be  related  to  the  EMP  plane  wave 
shielding  effectiveness  by  simple  correction  factors. 
(See  Chapter  2.) 


31 Document  Review.  Standard  Shielding  Effectiveness  Test 
Methods.  Technical  Memo  TM-19,  The  Boeing  Company, 
Huntsville,  AL  (U.S.  Army  Corps  of  Engineers,  31  May  1972); 
Measurement  of  Electromagnetic  Shielding  Effectiveness- A 
Review  of  Present-Day  Technology,  IRT  8194-017,  prepared 
by  IRT  Corporation  for  Defense  Nuclear  Agency  (16  January 
1980);  Recommended  Test  Procedure  for  the  Measurement  of 
Electromagnetic  Shielding  Effectiveness  of  High  Performance 
Shielded  Enclosures,  IRT  8494-016-1  (IRT  Corporation, 
21  March  1980);  R.  G.  McCormack,  Selection  of  Recommended 
Electromagnetic  /Radio  Frequency  Interference  Shielding  Effec¬ 
tiveness  Test  Procedures  for  Military  Tactical  Shelters,  ESL- 
TR-80-01 ,  prepared  by  the  U.S.  Army  Construction  Engineering 
Research  Laboratory  (Engineering  and  Services  Laboratory, 
Air  Force  Engineering  and  Services  Center,  January  1980). 


Using  the  appropriate  correction  factors,  the  approxi¬ 
mated  CW  shielding  response  curve  for  a  shielded  vol¬ 
ume,  and  Fourier  analysis,  the  shielding  for  plane  wave 
transients,  such  as  those  associated  with  HEMP,  can  be 
approximated. 

it  is  unlikely  that  any  radically  new,  low-cost  methods 
for  performing  simple  EMP  tests  on  shielded  structures 
will  be  developed  without  major  effort. 

In  view  of  these  points,  CW  testing  continues  to  be 
the  basic  low-cost  approach  for  assuring  continuing 
EMP  hardness  of  small  shielded  structures  such  as  shel¬ 
ters.  If  EMP  hardness  of  the  shelter  (with  all  internal 
equipment)  is  ascertained  after  it  is  manufactured 
and  assembled,  the  CW  test  approach  is  especially 
useful  for  measuring  any  shielding  deterioration.  Thus, 
having  concluded  that  CW  testing  will  remain  as  a 
compromise  method  of  testing  EMP  hardness,  the  test 
procedures  to  be  used  must  be  recommended.  In  a 
previous  study,  CERL  investigated  EMI  test  methods 
for  tactical  shelters.3*  The  study  concluded  that  IEEE 
299  test  methods  are  generally  preferable  to  those 
outlined  in  MIL-STD-285.  In  addition,  the  following 
changes  to  IEEE  299  were  recommended. 

1.  Where  possible,  use  electrostatically  shielded 
loops  for  the  small  loop  magnetic  field  tests.  If  these 
loops  give  an  inadequate  dynamic  measurement  range, 
unshielded  loops  may  be  used.  To  increase  measurement 
range,  the  unshielded  loops  should  be  used  with  mul¬ 
tiple  turns,  resonance  tuning,  and  impedance  matching. 

2.  For  general  use,  the  coaxial  loop  antenna  orienta¬ 
tion  is  preferable  to  the  coplanar  configuration  because 
it  can  detect  hidden  flaws  better  and  has  a  substantially 
greater  measurement  range.  The  coplanar  configuration 
is  recommended  if  its  measurement  range  is  adequate 
and  if  flaws  (other  than  seams)  can  be  visually  detected 
and  appropriately  tested. 

3.  The  number  of  small-loop  transmitting  antenna 
locations  should  be  increased  so  that  successive  test 
points  along  seams  are  never  more  than  2  ft  apart. 
(Test  points  around  penetrations  such  as  air  vents, 
power  line  filter  boxes,  and  signal  line  filter  boxes 


J*R.  G.  McCormack,  Selection  of  Recommended  Electro¬ 
magnetic  I  Radio  Frequency  Interference  Shielding  Effectiveness 
Test  Procedures  for  Military  Tactical  Shelters,  ESL-TR-80-01. 
prepared  by  the  U.S.  Army  Construction  Engineering  Research 
Laboratory  (Engineering  and  Services  Laboratory,  Air  Force 
Engineering  and  Services  Center,  January  1980). 
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should  not  be  more  than  1  ft  apart.)  If  this  becomes 
too  time  consuming  for  large  shelters,  the  large-loop 
test  of  IEEE  299  (or  a  commercially  available  seam 
leak  detector)  may  be  used  to  locate  the  points  of 
maximum  signal  leakage.  The  small-loop  test  may  then 
be  performed  opposite  these  points  to  determine 
worst-case  shileding  effectiveness. 

4.  The  small-loop  magnetic  field  test  need  not  be 
performed  at  frequencies  below  200  kHz  unless  the  shel¬ 
ter  has  specific  shielding  needs  below  that  frequency. 

5.  Selection  of  the  number  of  frequencies  at  which 
small-loop  magnetic  field  testing  is  to  be  done  in  the 
200-kHz  to  20-MHz  range  should  be  based  on  shielding 
criticality,  advance  knowledge  of  the  shelter's  future 
subjection  to  specific  electromagnetic  interference 
(EMI)  threats,  and  the  desired  overall  test  reliability. 

6.  Periodic  maintenance  testing  for  magnetic  field 
shielding  effectiveness  should  use  either  the  large-loop 
test  specified  in  IEEE  299  or  a  commercially  available 
seam  leak  detector  to  first  locate  points  of  maximum 
leakage.  These  points  should  then  be  tested  using  the 
small-loop  test  at  200  kHz  to  determine  worst-case 
shielding  effectiveness. 

7.  If  both  antennas  can  be  kept  well  away  from  other 
interfering  structures  when  the  ultra-high  frequency 
(UHF)  tests  (300  MHz  to  1  GHz)  are  performed,  the 
reference  signal  should  be  taken  using  a  calculated 
spacing  of  the  receiving  antenna  from  the  reflecting  wall. 
The  calculated  spacing  is  one-twelfth  of  a  wavelength. 

8.  Selection  of  test  frequencies  within  the  30-MHz 
to  10-GHz  range  should  be  based  on  the  desired  test 
reliability  and  on  engineering  judgment.  If  shielding  is 
not  highly  critical  to  the  shelter  mission,  a  single  fre¬ 
quency  within  this  range  is  sufficient;  a  test  frequency 
of  3  GHz  is  recommended  for  periodic  maintenance 
tests.  If  shielding  is  critical,  then  the  test  frequencies 
specified  in  MIL-S-55286C  (EL)  for  the  S-280  shelter 
are  recommended  (400  MHz,  1  GHz,  and  10  GHz). 
Furthermore,  if  shielding  is  critical  and  if  it  is  known 
that  the  shelter  will  be  subjected  to  high-level  EMI/RFI 
at  specific  frequencies,  then  additional  testing  should 
be  done  at  these  particular  frequencies. 

9.  A  horn  transmitting  antenna  may  be  used  for 
microwave  tests  if  it  is  needed  to  increase  the  measure¬ 
ment  range. 


Considering  the  findings  presented  in  Chapters  3  and 
4  and  CERL’s  experience  in  EMP  testing,  the  speci¬ 
fications  listed  in  IEEE  299,  with  the  modifications 
listed  above,  are  recommended  for  low-cost  methods 
of  assuring  continuing  EMP  hardness.  However,  it  is 
further  recommended  that  maximum  use  be  made  of 
either  the  IEEE  299  large-loop  test  or  seam  leak  testers. 
In  addition,  a  series  of  test  programs  should  be  pursued 
in  which  shelters  are  initially  submitted  to  threat-level 
EMP  tests  and  to  the  IEEE  299  tests  recommended 
here,  but  then  followed  by  extensive  seam  leak  or 
IEEE  299  large-loop  tests;  this  method  will  produce 
test  data  histories  showing  the  relative  merit  of  these 
low-frequency  approaches.  It  is  believed  that  if  EMP 
testing  is  used  as  a  base  with  corresponding  seam  leak 
test  data,  deterioration  in  EMP  hardness  will  be  demon¬ 
strated  by  a  deterioration  in  seam  leak  data. 


7  CONCLUSIONS 


This  report  has  evaluated  the  analytical  techniques 
applied  in  the  past  to  convert  CW  shielding  effectiveness 
test  data  to  values  meaningful  for  EMP  radiation.  Exten¬ 
sion  of  the  Shelkunoff  transmission  line  theory  of 
shielding  to  slotted  shielding,  as  proposed  by  Monroe, 
yields  inaccurate  results,  especially  for  test  loops.  A 
better  technique  is  the  Dipole-Moment  Approximation, 
which  gives  accurate  predictions  for  electrically  small 
slots.  Use  of  the  plane  wave  spectrum  approach,  as 
proposed  by  Villaseca,  along  with  Huygen's  principle, 
produces  results  that  are  in  good  agreement  with  the 
experimental  results. 

This  evaluation  has  also  produced  recommendations 
for  state-of-the-art  CW  test  techniques  for  evaluating 
the  EMP  hardness  of  tactical  shelters.  These  include: 

1 .  Use  the  specifications  listed  in  IEEE  299,  with 
modifications  as  listed  on  pp  30  through  31 . 

2.  Make  maximum  use  of  the  IEEE  299  large-loop 
test  or  seam  leak  testers. 
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APPENDIX  A: 

ELECTROMAGNETIC  COUPLING 
THROUGH  MULTILAYERED  SHIELDS  FOR 
AN  ARBITRARY  ILLUMINATING  SOURCE 


Introduction 

The  EMP/EMI  testing  of  shelters  requires  measuring 
the  coupling  between  two  test  antennas.  One  is  located 
external  to  the  enclosure,  and  the  other  is  placed  in  the 
interior  region.  Typically,  the  test  antennas  are  small 
loops  or  dipoles,  and  they  can  often  be  modeled  as 
elemental  magnetic  or  electric  dipoles.  While  formulas 


for  computing  the  coupling  between  two  coaxial  loops 
in  the  presence  of  planar  multi-layered  shields  can  be 
found  in  the  literature,29  general  formulas  for  aibitrary 
orientation  of  the  loop  or  dipole  antennas  are  not  readily 
available.  The  approach  documented  here  is  based  on 
the  plane  wave  spectral  representation;  this  method 
is  useful  for  calculating  the  shielding  effectiveness  of 
multi-layered  sheets  under  very  general  test  conditions. 
Since  the  method  is  capable  of  handling  arbitrary 
incident  fields,  loops  or  dipoles  oriented  in  arbitrary 
directions  can  be  accommodated. 

The  method  proceeds  by  considering  an  interface 
between  two  media  (see  Figure  Al)  and  expresses  the 
incident  fields  in  terms  of  a  plane  wave  spectrum.  The 
next  step  relates  the  fields  in  the  input  and  output 
terminal  planes  (arbitrarily  defined  on  two  sides  of  the 
junction  of  two  dissimilar  media)  via  a  transfer  matrix, 
T,  which  is  independent  of  the  nature  of  the  incident 
field.  The  elements  of  T  are  expressible  in  terms  of  the 
constitutive  parameters  of  the  two  media  and  the  loca¬ 
tion  of  the  terminal  planes.  The  multiple-junction  case 
is  handled  rather  simply  by  cascading  the  T-matrices 
for  the  individual  interfaces;  no  limitation  is  placed  on 
the  number  of  layers  in  the  shield. 


19R.  L.  Monroe.  EMP  Shielding  Effectiveness  and  MIL- 
STD-285.  HDL-TR-16 36/ AD77 1 997  (Harry  Diamond  Labora¬ 
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Validity  of  Applying  MIL-STD-285  to  EMP  Shielding  Effec¬ 
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System  Division  (Defense  Nuclear  Agency,  IS  April  1977); 
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magentic  Compatibility.  Vol  EMC-11,  No.  2  (May  1969). 
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OUTPUT  REFERENCE 

PLANE 

(PLANE  2) 


Figure  Al.  Media  interface  with  input  and  output  reference  planes. 
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The  following  section  presents  the  derivation  of  the 
T  matrix  that  relates  the  transformation  of  fields 
between  the  input  and  output  terminal  planes  1  and 
2  in  Figure  Al . 

Derivation  of  the  Transformation  Matrix,  T 

The  Cartesian  coordinate  system,  with  the  z-axis 
normal  to  the  interface, appears  to  be  the  obvious  choice 
for  expressing  the  fields  for  the  interface  problem  shown 
in  Figure  Al.  With  this  choice  for  axes,  there  are  four 
components  of  the  E-field  (EJ,  E$,  E$,  E^)  at  plane  1 
and  similar  ones  at  plane  2;  the  superscripts  a  and  b  refer 
to  incoming  and  outgoing  fields  that  must  be  related 
by  a  4  x  4  matrix.  However,  considerable  simplification 
results  if  an  alternative  choice  is  made  for  the  coordin¬ 
ate  system  in  which  the  various  fields  are  represented. 
Let  us  consider  an  incident  plane  wave  and  define  the 
plane  of  incidence  as  that  containing  the  wave  normal 
and  the  normal  to  the  interface.  The  planes  of  reflec¬ 
tion  and  transmission  can  be  defined  similarly.  We  then 
resolve  the  incident  E  and  H  fields  into  two  components: 
one  is  parallel  to  the  plane  of  incidence,  and  the  other 
is  perpendicular.  These  components  are  labeled  with 
subscripts  ||  and  1,  respectively.  Then  we  can  show  that 
these  two  components  propagate  independently  and 
remain  uncoupled  throughout  the  process  of  reflection 
and  transmission.  It  is  important  to  observe  that  if  we 
were  to  work  with  the  conventional  x,  y,  z  system,  there 
would  generally  be  coupling  between  the  various  com¬ 
ponents.  Thus,  for  convenience,  it  is  useful  to  introduce 
the  notation  where  the  superscript  a  indicates  incoming 
waves  into  the  planes  1  and  2,  while  the  superscript  b 
indicates  outgoing  waves.  Finally,  the  superscripts  1  and 
2  would  refer  to  the  fields  at  the  corresponding  ter¬ 
minal  planes  and  the  '  on  top  would  indicate  Fourier 
transforms. 

The  plane  wave  spectrum  representations  for  the 
various  field  components  take  the  form: 


Eu  =  //Eu  (a,0)ejla5<  +  *y-7'<l  +  d'>l  dad0  [Eq  Al] 
1  1 

E'b  =  jygib  (a  0  e ilex  +  -  v,(r  ♦  d,)J  da  d(j  [Eq  A21 

1  i 

Ej*  =  //Em  (a,  0)  e  +  dy  *  7,(r * dad|3  {Eq  A3} 
i  i 

E***  =  //Elb  (a,0)ej,ax+dy-7>(*  +  d»)l  dad0[Eq  A4] 
1  1 


where  yly 3  =  Vk^,  -  a5  -  0J,  k*  2  =  cuVi.i-  »nd  di>  da 


are  distances  of  the  terminal  planes  from  the  interface 
of  the  two  media. 

Similar  expressions  can  be  written  for  the  field’s 
parallel  components.  However,  these  fields  can  be  ana¬ 
lyzed  totally  independently  because  they  are  uncoupled. 

The  next  step  applies  the  curl  equations  to  obtain 
the  corresponding  magnetic  field  components  and 
applies  the  continuity  conditions  of  the  tangential 
fields  at  the  interface.  Instead  of  having  to  work  in  the 
spatial  domain,  one  can  apply  the  continuity  conditions 
equally  well  in  the  spectral  domain,  obtaining: 


[Eq  A5] 

E1*  ei7'd>  +  E,b  e~i7'd*  =  E“  ej7jd»  +  E*  e*j7ld» 

Vefc  [Eu  ej7'd>  -  E,b  e-j7>d>]  cosfl ,  = 

[E“  ej7>d» .  E*  e-j7»d’]  cos02 


where 


0,  =  angle  between  k,  and  the  normal  to  the  surface, 
vector  z, 


0,  =  angle  between  ka  and  the  normal  to  the  surface. 

k|-£  „  k,*z 

vector  z,  and  cos  0 ,  =  — -  cos  02  =  — — . 

Ikil  |k2| 


We  have  dropped  the  subscript  1  and  the  "  on  top, 
which  are  implicit. 

The  transfer  matrix  between  the  terminal  planes  1 
and  2  is  defined  as: 


{E1}  *  [T]  {E* }  [Eq  A6] 

i  1 

Following  the  steps  given  above,  one  can  obtain  the 
expression  for  the  transfer  matrix,  which,  for  the  1 
fields,  the  elements  of  which  are  given  by: 


.y^cotfi-v/ecotf  ,ejCT|d|_TA) 
•*■11  2Pj 


V^cos02  +  V<  COS0, 

2P, 


+  r.d.) 


T|  .>/£ cg«» ♦>/?«*.  jc, 

2P.  C 
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t,  *  V^cosOa-VicosO, 

A»  2Pj 

where  P2  =  %/ej  cos#  2. 

Similarly,  one  can  derive  the  transfer  matrix, 
T||,  for  the  parallel  component.  Its  elements  are 
given  by: 


TlI  ei(v,d,-r,<i,) 

lin  2Pj 

T(|  -  V^COSfl2-t-V^COSd,  »  7.d.) 

lln  2Pj 

T  -  Vei  cosfl2  +  y/ej  cost?  i  »Ti-i  +  7jdi) 


=  y5  cfg.-vgi6”*!  ei(,A  -  7,d,) 

Ha  2Pj 

For  the  multi-layered  case,  the  individual  T-matrices 
for  each  of  the  interfaces  may  be  cascaded  to  derive  a 
composite  T-matrix  for  the  entire  shield.  Then,  the 
transmitted  fields  in  the  extreme  right  region  can  be 
related  to  those  in  the  extreme  left  using  the  composite 
T-matrix.  Finally,  the  incident  fields  in  the  extreme  left 
region  can  be  assumed  known,  and  one  can  use  the  fact 
that  the  incident  field  is  identically  zero  in  the  extreme 
right  region.  Hence,  the  reflected  and  transmitted  fields 
in  the  extreme  left  and  extreme  right  regions,  respec¬ 
tively,  can  be  readily  calculated  from  the  knowledge  of 
the  composite  T-matrix  and  the  incident  fields. 

Determination  of  Incident  Fields 

This  section  briefly  describes  the  method  for  deter¬ 
mining  the  fields  generated  by  elemental  electric  and 
magnetic  dipole  sources.  We  consider  four  different 


1 .  Electric  dipole  oriented  along  the  z-direction 

2.  Electric  dipole  oriented  along  the  x-direction 

3.  Magnetic  dipole  oriented  along  the  z-direction 


Electric  Dipole  Oriented  Along  Z-DIrection 

For  a  z-oriented  dipole  source,  the  z-component 
of  the  electric  vector  potential  At  satisfies  the  wave 
equation: 

7JA,  +  k’A,  =  -Idl  5(x-x')  6(y-y')  6(x-x')  [Eq  A7J 

where  Idl  is  the  dipole  moment  of  the  elemental  dipole. 
This  equation  for  At  can  be  solved  using  the  transform 
technique  to  yield: 


A*(x,y,z)  =  2(2wyjfm 


.  e  il«(x-*V  ff(y-y’)  -  s/k1  -o’  -ff‘  li-il'ldodfl 

From  Maxwell’s  equations,  we  have: 

dA»  +1  d*A» 

E„  = : —  [Eq  A8] 

3y  ^  ](oe  dzdx 

9  A,  1  3*AX 

The  expressions  for  the  various  components  of  electric 
and  magnetic  fields  are  given  in  Eqs  A10-A14. 

[Eq  A10] 


[Eq  All] 


Hy  =  2(2n?IL^^e>[  ,dad0;U 


[a(x-x’)+0(y-y')-7l*-z,l] 


1  <to‘V  |Eq  A13] 


4.  Magnetic  dipole  oriented  along  the  x-direction.  c  —  /T*  .»  |  m  re.  iui 

*  (joie)2(2pi)1Ji-VkJ^i  *  ^  1 


All  of  these  problems  are  most  conveniently  solved  in 
the  transform  or  spectral  domain,  where  the  wave 
equations  for  electric  or  magnetic  vector  potentials 
become  algebraic. 


Electric  Dipole  Oriented  Along  x-Direction 

For  a  dipole  oriented  along  the  x  direction,  the  field 
components  are  given  by: 
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(Eq  A15J 


Hy  -  ejl“(x-x’>+«sr-»’>-1'|l-*'l»da«J|3 

H*  =  Wf-  eJM  ^  ^  tEq  A161 

n^Aiu 


^  *  (jue)2(2n)1Jj-m  ejM  ^  ^  {EqA18] 


El  =  (j^^/;-:jttei'  '  dad^ 


[Eq  A19] 


Magnetic  Dipole  Oriented  Along  the  z-Direction 

For  a  z-oriented  magnetic  dipole  of  moment  Kdl, 
we  solve  first  for  the  magnetic  vector  potential  Fz  and 
derive  the  E-  and  H-fields  from  this  potential.  The  Held 
expressions  are: 

[Eq  A20J 

[E,A2“ 

(Eq  A22] 


H'  “  (jun)2(2n)‘ 


(EqA24J 


Magnetic  Dipole  Oriented  Along  the  x-Direction 

Finally,  for  an  x-oriented  magnetic  dipole  of  moment 
Kdl,  the  corresponding  expressions  are: 


[Eq  A25J 


Ey  -  /JLI ei,a(*-*')+tf(,-y’)-7l,-,'lldo  dff 


Ez  =  2(2^  ^  eH  lda<^  [Eq  A26) 


«■  ’  ss§k?£  •" 1 "» •*  » 


H”  “  (wna,)1!-  Vk5?  *11 '  ‘k"dii  16,1  A2S| 


Once  again,  one  can  work  initially  in  the  transform  do¬ 
main  to  carry  out  the  multiplication  of  the  T-matrices 
and  use  the  incident  field  components  to  derive  the 
transmitted  fields,  also  in  the  transform  domain.  These 
can  then  be  Fourier-inverted,  and  the  field  coupling  into 
another  dipole  at  the  receiving  end  may  be  computed 
by  taking  a  dot  product  of  the  appropriate  field  com¬ 
ponent  with  the  electric  or  magnetic  dipole  moment. 

Coordinate  Transformation 

Let  us  consider  the  problem  of  coordinate  transfor¬ 
mation  from  the  x,  y,  z  coordinate  system  to  the  system 
defined  by  the  orthogonal  triad  consisting  of  the  wave 
normal  and  the  vectors  parallel  and  perpendicular  to 
the  plane  of  incidence.  To  accomplish  the  transforma¬ 
tion,  we  use  a  rotation  matrix  R.  For  the  E  fields,  the 
matrix  R  is  defined  by: 


[Eq  A30) 


— 

-sind 

cosd 

0 

[Eq  A23] 

R  = 

-cosd  cos# 

-cosd  sin# 

sind 

[Eq  A24] 

sind  cos# 

sind  sin# 

cosd 

[Eq  A31] 


where: 


V  W 

k-S 

/|kia  -  (k*£)J 
k-x 

/|k|2  -  (k-z)3 
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APPENDIX  B: 

COMPUTER  PROGRAMS:  MULSH,  APRAD, 
SLOT,  SLITCP,  and  SLITCA 


PROGRAM  SL0T2 l INPUT .OUTPUT, TAPE5- INPUT ,TAPE6=QUTPUT) 

COMPLEX  RKD,RKL,RKDPO,RKLPO 

COMPLEX  ZSL,rOP,2tl,ZDT,ZDB,Zl,2Lf  ,ZLB 

Pl-3. 1  415926536 

URITE  l&, 100 ) 

100  FORMATS, 2X,*LARGE  SLOT =01 , SMALL  SL0T=00*> 

READ ( 5 , 200 ) IS1 Z 

200  F0RHAT  < 12 ) 

IF ( IS IZ  .EQ.  1)  GO  TO  10 
SL=0.25 
A-7.9375E-4 
GO  TO  20 
10  SL=0.5 

A=1 .5875E-3 
20  CONTINUE 
R=.3048 
WRITE  <6 , 101 > 

101  F0RNAT(/,2X,*ENTER  START  FREQ.  (MHZ)*) 

READ (5,201 )FSTART 

201  FORMAT (F  14.2) 

FR=F3TART 
WRITE (6 , 1 0?) 

107  FORMAT (/ ,2X,*ENTER  MAX  FREQ.  (MHZ)*) 

READ(5,203)FRMAX 
203  FORMAT  <F 1 4.2) 

FRINC-10.** ( . 125) 

IFLP=0 

1234  F=FR*1.E+4 
1FLP= IFLP+ 1 
ALAN=3.0E08/F 
BETA~2. *PI/ALAM 
TBL=2.*BETA*SL 
FBL=4.*BETA*SL 
SOTBL'O .0 
S0FBL=0.0 
SITBL~Q.O 
SIFBL=0.0 

C  CALCULATE  SQTBL 
M-0 

30  N'N+ I 
I7N=2*N 

CALL  FACT(ITN,ITNF) 

501  FORMAT (14,14) 

TERM-(-1 . )*( (-1 . )**N)*( (TBl )**ITN)/ ( ITH*ITNF) 
SOTBL=SOTBL+TERM 
IF <N  .LT.  5)  GO  TO  30 
C  CALCULATE  SOFBL 
N=0 

40  N-N+1 
ITN--2*N 

CALL  FACHITN.ITNF) 

TERH=(-1 .  )*t (-1 . )**N)*( <FBL)**ITN)/(ITN*ITNF) 

SOFBL-SOFBL+TERM 

IF i N  .LT.  5)  GO  TO  40 


t.  CAllULAIt  :ii  fbl. 

N  ■-  I 

50  N-H-t  I 

ITNPQ=2+N+1 

CALL  FACT ( ITNP0,1TNPQF) 

TER«=  ( ( -1 .  (TBL )* *ITNPO »/ < ITHPO=*ITNPOF ) 

S1TPL=SITPL+TER« 

IF(N  .LT.  4)  GO  TO  50 
C  CALCULATE  SIFBL 
N=-f 

60  N=N+1 

1TNP0=2*N+1 

CALL  FACTUTNPO,  ITNPQF) 

TERN=(  (-1 .  ):**N)*(  (FBL)**ITNPO)/( ITNPO*ITNPQF ) 
S I FBL =S I F  BL  +  TERM 
IF(N  .LT.  4)  GO  TO  60 
RAB=15.*( (2.+2.*C0SUBL) )*SOTBL 
4-COS  <  TBL  >*S0FBL-2.=*S3M<  TBL  )*SITBL 
J+SIN1  TBL ) *S I FBL > 

20=120. *<AL0G<SL/A>-1 .0-.5*AL0G(2.*SL/ALAN) ) 

GAM=2.0*RAB/Z0 

EPG=EXP<GAM) 

ENG=1./EPG 

EPTG=EXP(2.*GAN) 

Ei1TG=1  ./EPTG 
SH7G=<EPTG-ENTG)/2. 

CHG=(EPG+EM6)/2. 

CHGS-CHG**2 

C0SSBL=(C0S(FE7A*SL>)**2 
502  F0RHAT(E14.7,E14.7) 

RCD=  <  ZO/ 2 . ) *  <  SHTG/ ( CHGS-COSSBL ) ) 

XC0=(Z0/2 .  )* ( (-1 . )*SIN(TBL)/< CHGS-COSSBL) ) 

RCDS=RCD**2 

XCBS=XCD*»2 

TOP=CMPLX (RCB ,  ( -1  - )  •* X C D ) 

BOT=RCDS+XCDS 

ZSL=<<<3?7.)**2)/4.)*<T0P/B0T> 

ZMAG=CABS ( ZSL ) 

102  FORMAT (/, 2X,*ZSL=* ,E1 4.4 1 
C  CALCULATE  BIPOLE  IMPEDANCE 
ZRDT=1.-<BETA**2)*<R**2> 

ZIDT*BETA»R 

ZRDBM-1.)*(BETA**2)*(R**2) 

Z I  llB  ~  Z  IDT 

ZDT=CMPLX<ZRDT,ZIDT) 

ZBB=CMPLX(ZRDB,ZIDB> 

ZD=(377.)*<ZDT/ZDB) 

ZDNAG=CABS(ZD) 

C  CALCULATE  LOOP  IMPEDANCE 
ZRLT=ZRDB 
ZILT=ZIDT 
ZRLB=ZRDT 
ZILB=ZIDT 

ZLT=CMPLX(ZRLT,ZILT) 

ZLB=CMPLX(ZRLB,2ILB) 

ZL=<377.)*<ZLT/ZL8) 

ZLMAG=CABS(ZL) 
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104  FORMAT (/ , 2X, IE1 4 .4) 

C  CALCULATE  SHIELDING  EFFECTIVENESS  FOR  DIPOLE 
RKD=2D/2SL 
RKDPO=RKBH.O 
DKN=CABS(RKB) 

DKNO'CABS(RKDPQ) 

S£D=20. *AL0610<BKNO**2/(4.*DKN) ) 

C  CALCULATE  SHIELDING  EFFECTIVENESS  FOR  LOOP 
RKL=ZL/ZSL 
RKLP0=RKL+1  .0 
EKM=CABS(RKL) 

EKNO-CABS(RKLPQ) 

SEl=20.*AL0G10<EKHQ**2/<4.*EKM) ) 

IFlIFLP  .GT.  1)  GO  TO  7? 

URITE(6,105> 

105  FORMAT (20X,*FREQ(HZ)*, 10X, *SE  DIPOLE*, 10X,*SE  LOOP*,//) 
77  CONTINUE 

UNITE (6, 106)F,SED,SEL 

106  FORNAT (18X,E10.3,8X,E10.3,9X,E10.3> 

FR=FR*FRINC 

IF1FR  .LE.  FRHAX)  GO  TO  1234 
URITE(6,103) 

10J  FORMAT!/, 2X,»C0NTINUE?  T=I,  N=0*> 

READ(5,202)IROT 
202  FORMAT (12) 

IF ( 1R0T  .EO.  1)  GO  TO  20 

STOP 

END 

SUBROUTINE  FACT (N.NFACT) 

NFACT=N 

IF(N  .EQ.  0)NFACT=1 
H=N 

10  continue: 

M=M-1 

IF (M.EO. 1 ) GO  TO  20 
IF (M .EU .0 )G0  TO  20 
NFACT=NFACT*M 
GO  TO  10 
20  CONTINUE 
RETURN 
END 
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PROGRAM  SL I  TCP ( INPUT, OUTPUT, TAPE5= INPUT, TAPE6=GUTPUT ) 
REAL  LAM,HR£FR,HR£FI ,HREFM,HRI,HRR,HTR,HTI ,HDTZM 
COMPLEX  HR,HTHE,HX, ARG,HDZ,HDZT 
DIMENSION  XP(5),DXP<5) 

PI=3.1415?26536 

UR  I T E  <  6 , 102) 

102  FORMAT*/, 2X,*LARGE  SL0T=01 ,  SMALL  SLOT=00*) 

READ<  5 , 202 ) 1  SLOT 

202  FORMAT (12) 

WRITE (6,100) 

100  FORMAT </,2X,*ENT£R  START  FREQ  (MHZ)*) 

R£AD<5,200)FR 

200  FORMAT(EIO.I) 

UR  I TE (6,101 ) 

101  FORMAT (/,2X,*ENT£R  MAX  FREQ  (MHZ)*) 

READ(5,201 )FRMAX 

201  FORMAT (E10. 1 ) 

CELT-. 125 
FRINC=10.**DELT 
UR ITE ( 6 , 1 03  > 

103  FORMAT(//,30X,*FREQ  (HZ)*, 10X,*SE  (DB)*,//) 

C  BEGIN  FREQ  LOOP 

10  F=FR*1 . E+6 
LAM=3.0E08/F 

C  CALCULATE  REF  H  FIELD  ON  AXIS  AT  3  FEET 
BETA=2.*PI/LAM 
R=.91 AA 
BETAR=BETA*R 

HREFR=C0S(BETAR)-SIN(BETAR)/BETAR-C0S(BETAR)/(BETAR**2) 
HREFI=SIN(BETAR)/(BETAR**2)-SIN(BETAR)-C0S(BETAR)/BETAR 
HREFR’HREFR*(8ETA**2)/( A. *PI*R) 
HREFI=HREFI*(BETA**2)/(A.*PI*R) 
HREFM=SQRT(HREFR**2+HREFI**2) 

IF ( ISLOT  .EQ.  1)  GO  TO  20 
IFUSLOT  .EO.  0)  GO  TO  30 
20  XP  < 1 )  = . 25 
XP<2)=.5 
XP(3)=-.25 
XP(A)=-.5 
XP(5)=0.0 

BXP< 1 )=BXP(3)=DXP(5)=.25 
DXP(2)=DXP(A)=.125 
DY=. 0015875 
EL=1 .0 
AMX=. 02132 
60  TO  AO 
30  XP(1 )=.125 
XP(2)=.25 
XP<3)=-. 125 
XP(A)=-.25 
XP(5)=0.0 

DXP(1)=DXP(3)*DXP(5)=.125 

DXP(2)'DXPi,A)*.0625 

DY~. 00079375 

EL=0.5 

AMX=. 002665 
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40  CONTINUE 

HD2T= (0. 0,0.0) 

DO  2000  1=1,5 

C  CALCULATE  HX  ALONG  SLOT 
Z=.45?2 
X=XP(I) 

DX=DXP<1) 

ALPHA=ATAN(X/.4572) 

THETA=PI/2.-ATAN (X/.4572 ) 

R=. 4572/COS (ALPHA) 

BETAR-BETA*R 

HRR=<C0S<THETA)/(2.*PI>)*(<C0S(BETAR)/<R*+3>) 
++6ETA*SIN(BETAR)/(R**2) ) 

HRI  =  ( COS(  THETA) /<2-  *PI >  > ♦  ( BET A*CQS < BETAR > / < R=r*2 ) 
+-SIN(BETAR)/(R**3)J 
HR=C«PLX(HRR,HRI) 

HTR=COS(BETAR)-SIN(BETAR)/BETAR-COS (BETAR)/ (BETAR+*2) 
HTR=HTRt(-1.0)*(BETA**2>*SIN(THETA)/<4.*PI*R) 
HTI=SIN(BETAR)/(BETAR**2)-SIN(BETAR)-CQS(BETAR)/BETAR 
HTI=HTI*(~I . 0)*(BETA**2)*SIN( THETA)/ (4.+PI*R) 
HTHE=CHPLX(HTR,HTI ) 

HX=(HR*X)/SaRT<X4*2-K.45?2)**2> 

♦+HTHE*(. 4572) /SORT <X**2+(. 4572)4*2) 

HX=HX*DX 

C  CALCULATE  CONTR.  TO  SHAOOU  SIDE  FIELD 
RHRP=SQRT(X**2+(.4572>**2) 

RHRPH=(-1 .0>*RMRP 
ARG=CHPLX(0.0,RHRPM> 

HB2=<S.<A«X*HX*((.4572)**2)*CEXP(ARG)/((RNRP**5)+4.*PI) 
2000  HDZT=HDZT+HDZ 

HDZT=HDZT/(5.0*EL) 

HDZTH~CABS(HDZT) 

SE=20.*ALQG10(HREFrt/HDZT«> 

URITE(6,104)F,SE 
104  F  ORHAT ( 30X ,E10.3,5X,F10.3) 

FR=FR*FRINC 

IF(FR  .LE.  FRHAX)  GO  TO  10 

STOP 

END 

C> 
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PROGRAM  SLITCA< INPUT,GUrPUTfTAPE5-INPUI , TAPEA=QUTPUT > 
REAL  L AM, HREFR, HREF  I, HREFH.HR I, HRR,H'fR,HTI fHDTZH 
COMPLEX  HR,HTHE,HX, ARG,HDZf HD2T 
DIMENSION  XP(4) ,DXK4 ) 

FI=3.14I592A536 
URITEU,  102) 

102  FORMAT (/ ,2X,*LARGE  SL0T=01,  SMALL  SLDT=00*) 
READ(5,202)ISL0T 

202  FORMAT (12) 

URITEIA, 100) 

100  FORMAT (/,2X,*ENTER  START  FREQ  (MHZ)*) 

READ(5,200)FR 

200  FORMAT <E 10.1 ) 

UR  I TE  <6.101 ) 

101  FORMAT </, 2X,*ENTER  MAX  FREQ  (MHZ)*) 

READ(5,201)FRNAX 

201  FORMATIEIO. 1 ) 

OELT= . 1 25 

FRINC=10.*:#DELT 

UR1TE(6.103> 

103  FORMAT  (II , 30X , *FREQ  <H2>», 10X,*SE  < DB )=«,// ) 

C  BEGIN  FREQ  LOOP 

iO  F=FR*1.E+6 
LAM-3. 0E08/F 

C  CALCULATE  REF  H  FIELD  ON  AXIS  AT  2  FEET 
BETA=2.*PI/LAM 
R=.A1 

BETAR=BETA*R 

HREFR=C0S(BETAR)/(R**3)+(BETA*SIN(BETAR) )/ (R**2) 
HREFI=<-1.)tSIN(BETAR)/(R**3>+<BETA*C0SlBETAR))/<R+*2) 
HREFR=HREFR/<  2.*PI ) 

HREFI=HREFI/(2.*PI) 

HREFM  -SORT (HREFR**2+HREF 1**2) 

IF ( ISLOT  .EQ.  1)  GO  TO  20 
IFdSLOT  .EQ.  0)  GO  TO  30 
20  XP( 1 )=.25 
XP(2)=.5 
XP(3)---.25 
XP ( 4 )~— . 5 
DXP ( 1 )=DXP(3) -.25 
DXP(2i=DXP(4>=.125 
DY~. 0015875 
EL=1 .0 
AMX=. 02132 
GO  TO  40 
30  XP<1)~.!25 
XP \ 2 ) = . 25 
XP(3)=-.125 
XP(4)--.25 
DXP ( 1  )=DXP(3>". 125 
DXP(2)=DXP(4 )=.0625 
DY= .00077375 
El-O. 5 
AMX=. 002665 
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40  CONTINUE 

HDZT“(0. 0,0.0) 

DO  2000  1=1,4 

C  CALCULATE  ttX  ALONG  SLOT 
Z  =  .305 

x=xpm 

DX=DXF ( I ) 

THETA=ATAN(X/.305) 

R=. 305/COS (THETA) 

BETAR=BETA*R 

HRR= (COS<  THETA)/<2. *PI) )*( (CQS(BETAR)/(R+»3) ) 
++B£TA*SIN(B£TAR)/(R**2)) 

HR I=( COS ( THETA) /<2.*PI) )*(BETA*C0S(BETAR)/(R**2) 
+-S1N(BETAR)/(R**3)) 

HR=CHPLX(HRR,HRI) 

HTR=CQS ( BETAR) -SIN (BETAR) /BE  TAR-COS < BET AR)/(BETAR**2> 
HTR=HTR+( -1.0)* (BETA**2)*SIN< THETA )/<4.*PI*R) 
HTI=SIN(BETAR)/(BETAR**2)-SIN(BETAR)-CPS(BETAR)/BETAR 
HTI=HTI*(-1 .0)*<  BETA**2 )  *SIN(  THE  TA)./(4.*PI*R) 
HTHE=CHPLX(HTR,HTI) 

HX=<HR+X)/SQRT(X**2+(.305)**2) 

+tHTHE*X*( .305) /SORT ( iX**2)*(X**2+( .305)**2) ) 

HX=HX*2.*DX 

C  CALCULATE  CONTR.  TO  SHADOU  SIDE  FIELD 
RHRP=SQRT(X**2+(.305)++2) 

RHRPH=i-1 .0  )*RHRP 
ARG=CNPLX(G.0,RrtRPf1) 

HDZ=6.  *ANX:»HX*(  - 1 .0)  *X* ( .305):*CEXP(ARG)/(  (RrtRP**5)+4. *PI ) 
2000  HDZT=HDZT+HDZ 

HDZT=HDZT/(5.0*EL) 

HDZTN=CABS(HDZ7) 

SE=20.*ALDG10(HREF«/HDZTN) 

UFiITE(6,104)FfSE 
104  FORNAT(30X,E1 0.3,5X,F10.3) 

FR=FR*FRINC 

IF (FR  .LE.  FRNAX;  GO  TO  10 
STOP 


PROGRAM  MULSH  (INPUT, OUTPUT, TAP£5=INPUT ,TAPE6=0UTPUT> 

DIMENSION  F REQ < 1 00) ,FLOG( 100) ,SHEF( 100) , BCD < 3 ) ,BF (5000 ) ,SH( 100, )0) 
DIMENSION  LX(7),LY(7),LM(7),LS(7) 

COMNGN/PARAMS/NOPK?) ,D(7) ,SIGR<?> ,EHUR(7) ,EPSR(7) ,DLUP(7) ,FRCN(7) 
1 ,NG0,LL(7) 

COMMON/CONSTS/PI ,EMU0,EPS0,SIGC 
COMMON/STPCOM/IMAGE (2100) 

COMPLEX  EJ,  SUND,  CEP,  CEN,  PROD,  T2MD,  CE,  FACT 
COMPLEX  ZIN 

COMPLEX  GS<7>,  TAU(7),  PSY<?>,  HT(7),  SZR(7),  P0LE(7),  CSPL(7) 
COMPLEX  TRAN(IOOO),  REFL(IOOO) 

C 

PI=3. 1415926536 
TP=2.*PI 
ERR=1.E-10 
EMU0=1 .26E-6 
EPSQ=8.854E-12 
SIGC=5 .80E+7 
VL=3.E+8 
EJ=CMPLX (0. , 1 . ) 

C 

C . BEGIN  PROBLEM  LOOP 

KK=0 

1000  CONTINUE 
KK-KK+1 

URITE (6 ,100)KK 
CALL  IADATA 
MC=NOPT ( 1 ) 

MCM=HC-1 

IHAX=N0PT(2) 

IF ( IMAX.GT . 1 000) STOP  1000 
IDBP=NQPT (3 ) 

IPLT=NQPT (4) 

C 

A  =DLUP< 1 > 

DAL  =  DLUP  <  2 ) 

DD  =  0. 

DO  5  N=2,NCN 
5  DD=DD+D(N) 

ZMD  =  D( 1  )+B(HC) 

RDR  »  ZMD+DD 
C 

C . NUMERATOR  INTEGRAL 

RS;A*A+RDR*RDR 
R=SQRT (RS) 

TERN=EMUR(1)*R*R*R/A 

C 

C . BESSEL  TERMS 

AL=0. 

DO  10  I =1,1 MAX 
X*AL*A 

CALL  BESJ(X,1 ,BJ,ERR,IERROR) 

IF (IERRGR.GE. 1 .AND. I.GT. 1)G0  TO  70 
BF(I)=BJ 
10  AL=AL+DAL 
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1 005  CONTINUE 

C . BEGIN  FREQUENCY  LOOP 

FR  =  FRCNil) 

FRHAX  =  FRCN(2) 

FRINC  =  10.*+FRCN<3) 

11=0 

1010  CONTINUE 
11=11+1 

C 

F=FR*1.E+4 

U=TP*F 

C 

C . RESULTANT  PROPAGATION  CONSTANT ,  GANNA,  FOR  EACH  LAYER  AT  THIS 

C  FREQUENCY . 

DO  12  H  =  1 ,MC 
ENU  =  £NURtH)+EHUQ 
EPS  =  EPSR<H)*EPSO 
SIG  -  SIGR(H)*SIGC 
GNR  =  -  U*U*EHU*EPS 
GNI  =  U+ENU+SIG 
12  GS<N)=CNPLX(GHR,GHI) 

C 

C . INTEGRATION  LOOP 

AL=0. 

SUHB=0. 

HT (NC)=Q. 

SZR(NC)  =  1. 

DO  50  1=1 , IHAX 
ALS=AL»AL 
DO  14  N  =  1 ,NC 

14  TAU(N)  =  CSQRT < ALS  +  GS ( H > > 

C 

C . HYPERBOLIC  FUNCTION  CALCS 

DO  18  H=2,HCH 
PSY(N)=TAU(Hi*D(N) 

CEP  =  CEXP  < PS Y ( H ) ) 

CEH  =  1 ./CEP 
CSPL(N)=(CEP+CEH)/2. 

HT  <H)  =  (CEP-CEH)/ (CEP+CEN) 

18  CONTINUE 

C 

C . CALCULATE  IHPEDENCE  RATIOS . 

DO  32  NR  =  1,  NCN 
N  =  NCN-NR  +  1 
HP  =  N+1 

32  SZR(H)  =  (TAU(N)/TAU(HP) )*(S2R(NP)+HT(NP) )/( 1 ,+SZR(NP)*HT(HP) ) 
1  *<EHUR(HP)/EHUR(H) ) 

C 

C . CALCULATE  ALL  POSSIBLE  POLES 

POLE ( 1 )  =  1.  ♦  SZR < 1 ) 

PROD  =  2 . /POLE ( 1 ) 

DO  34  H  =  2,  NCN 

POLE(N)  =  I.  ♦  S2R(H)*HT(H) 


34  PROD  *  PROD/(POL£(H)*CSPL(H)) 

TRAN(I)  =  PROD 
ZNAG=CAfiS<SZR(2> ) 

ZR£AL=REAL  <SZR(2) ) 

ZIH=AIHAG(SZR(2)> 

ZPHA=ATAN2 (ZIH,ZREAL ) 

500  F0R«AT(/,2X,*SPACIWG=*,F8.5f5X,*MAG=-SEH.7,5X,*PHASE=*,E14.?) 

C  REFL(I)  =  (SZR(1)-1.)/(SZR(1)+1.) 

TZHD  =  TAU( 1 )*ZMD 
CE=CEXP(-TZHB) 

FACT  =  ALS*TRAN(I)*CE*BF(I)/TAU(1 ) 

40  SUHD  =  SUHD  +  FACT 

C . DEBUG  PRINTOUT _ _ 

IF(1DBP. £0.0)60  TO  42 

IF(HOD(I,100).NE.O)GO  TO  42 

WRITE<6,114HGS<N),TAU<H),PSY<H>,HTiH>,H=1,HC> 

WRITE (6, 1 15) (CSPL (N) , SZR(H) ,POLE(H) ,N=1 ,NC ) 
WRITE(6,116)ALS,PR0D,CE,BF(I> 

WRITE(6,11?)FACT,SUHD 
42  CONTINUE 
AL-AL+DAL 
789  CONTINUE 

50  CONTINUE 

C  . END  OF  INTEGRATION  LOOP . 

C 

51  CON  i  INUE 
?'JP>SUND*DAL 
CSUHD-CABS(SUMD) 

SHEF(II)=99?9. 

IF(CSUND.GE.1.E-40)SHEF(II)— 20.*ALOG1O(TERN*CSUHD) 

SH ( 1 1  ,KK ) =SHEF ( II ) 

FREQ ( I I ) =F 
FLQG(II)~AL0G10<FR) 

FR=FR*FRINC 

IF(FR.LE.FRHAX)GO  TO  1010 
C . END  OF  FREQUENCY  LOOP 

C 

URITE(6,104) 

WRITE (6, 105) (FREQ  < I ) ,SHEF  ( I ) ,1=1,11) 

WRITE(6, 1 06) 

IF (NGO.EQ. 1 ) GO  TO  1000 
C 

C . PRINTER  PLOT 

IF(IPLT.LE.0)60  TO  66 
BCD ( 1 )=1 OHSH I ELDING 
BCD(2)-10HEFFECTIVEN 
BCD(3)=10HESS 

CALL  STPLT2 ( 1 ,FLOG,SHEF, II ) 

CALL  PLQT1 (0,2,20,8, 10) 

CALL  PL0T2(INAGE,2. ,-5. ,200.,0. ) 

CALL  PLOT 3 ( 1 H . ,FLOG,SHEF,II, IHAGE) 

CALL  PL0T4(28, BCD, IHAGE) 

WRITE<6,107) 
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£ . END  OF  PROBLEM  LOOP 

C 

c 

C . DRUM  PLOT 

CALL  CONPRS 
CALL  BGNPL(O) 

DO  65  K=1,KK 
DO  61  1=1,11 

61  SHEF(I)=SH(I,K) 

IF (K.GT.1 ) GO  TO  62 
SHI =1 00 . 

SH2= 1 .3 
SV1 =0 . 

SV2=40. 

DO  60  1=1,7 

60  LX( I )=LY( I )=LH(I )=10H 
LX<3)=10HFR£QU£NCY 
LX(4)=1QH(H2) 

LY(2)=10HSHIELDING 
LY(3 )=10HEFFECTIVEN 
LY (4)=10HESS  < DB > 

CALL  TITl£(LL,-70,LX,70,LY.40,8.,5.) 

CALL  XL0G(SH1 , SH2 , S0 1 , S V2 ) 

62  CALL  CURVE ( FREQ ,SHEF, II ,12) 

65  CONTINUE 

CALL  ENDPL(KK) 

CALL  DONE PL 

66  CONTINUE 
C 

GO  TO  71 

70  UR I  IE (6, 1 12>  IERROR 
STOP  10 

71  CONTINUE 

100  FGRMATOHI ,/////,  30X , *PROBLEN  *,I2,*  INPUT  DATA*,//) 

101  FORMAT (20X,*FR£QUENCY  =  *,EI3.4,+  HZ*,1QX,*LAHDA  =  »,E13.6) 

102  FORMAT (  5(5X,2E10.3> > 

104  FORMAT! 1  HI ,30X,*FREQ(HZ >*,1 OX, *SHEF(DB )*,//) 

105  FORMAT (30X,E10.3,5X,F10.3) 

106  FORMAT ( 1 HI ) 

107  F0RNAT<30X, ♦FREQUENCY  (MHZ)  EXPONENT  *,1H1) 

109  F0RMAT(25X,2E13.6,5X,E13.6,5X,2E13.6) 

110  FORMAT!//) 

111  F0RNAT(//,20X,*DEN0NINAT0R  TOO  SMALL  FOR  F=*,E13.6,+  HZ*,//) 

112  FORMAT UOX,*BESSEL  FUNCTION  ERROR,  IER  =  +,I4) 

114  F0RMAT(4(4X,2E13.6) ) 

i  15  F0RMAT(3(2X,2EI3.6)) 

116  F ORMAT </,5X,E13.6,2(5X,2E13.6),5X,E13.6) 

117  FORMAT!/ ,25X,2( 1 0X,2E13.6) ) 

END 

SUBROUTINE  1ADATA 

COMMON/PARAMS/NOPH?) , T (?),SIGR(?) ,EHUR(7) ,EPSR(7 ) ,DLUP(7) ,FRCN(7) 
1 ,NG0,LL(7) 
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f 


5 


to 


501 

502 

503 

504 

505 

506 
601 
602 

603 

604 

605 


C 


NGG^O 

WRITE (6, 501 ) 

READ. 5,601 )N0PT(1) 
KONGf-mi 
T(1>=T<KC)-.305 
DO  5  IU=1,NC 
EPSR(IU)=1 . 

ENUR< I U )  =  1 . 
CONTINUE 
HOP  T  <  2 )  =  1 000 
N0PT<3)=0 
NOPT ( 4  >*0 
DLUP ( 1 )= . 152 


DLUP(2)=. 1 
SIGR(1)=SIGR(HC)=0. 

HCH=NC-1 
DO  10  IKI=2,MCN 
URITE(6,502)IKI  a 

READ(5,602)T (IK1)  f 

URI TE ( 6 ,503 )  IKI  • 

READ(5,603)SIGR(IKI) 

CONTINUE 
WRITE (6,504) 

READ(5,604)FRCN( 1 ) 

URITE(6,505) 

READ ( 5,605 )FRCN( 2) 

FRCN<3)=. 125  ^ 

URITE(6,506)  ** 

READ ( 5,601  )NGO  ! 

FORMATS, 2X,*ENTER  NO.  OF  LAYERS*) 
F0RMAT(/,2X,*ENTER  THICKNESS  OF  LAYER*, 2X, 12) 
FORMAT (/,2X,*EN TER  SIGR  OF  LAYER*, 2X, 12) 
FORMAT ( /,2X, *ENTER  START  FREQ.*) 

FORMAT (/,2X,*ENTER  MAX  FREQ.*) 
F0RMAT(/,2X,*M0RE  DATA  LATER?  Y=1,N=0*) 

FORMAT (12) 

FQRMAT(F1 4.8) 

FORMAT (F14.8) 

F0RMAT(F14.8) 

FORMAT (F14.8) 

RETURN 

END 

SUBROUTINE  BESJ(X,N,BJ,D,IER) 


BJ-.O 

IF(N) 10,20,20 
10  I ER= 1 


RETURN 

20  I F ( X )  30 ,30 ,31 
30  IER=2 


RETURN 

31  IF  <  X- 1 5 . )32,32,34 

32  NTEST=20.+10.*X-X**2/3 
GO  TO  36 
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34  NTEST-90.+X/2. 

36  IFtN-NTEST  >40,38,38 
38  IER-4 
RETURN 
40  IER=0 
N1=N+1 
BPREV=.0 

(, 

C  COHPUTF  STARTING  VALUE  OF  N 

C 

IF (X-5. ) 50, 60, 60 
50  NA=X+6. 

GO  10  70 

60  HA=1.4*X+60./X 
70  «B=N+IFIX(X)/442 
rtZERO=NAXO (rtA  ,NB ) 

C 

C  SET  UPPER  LIHIT  OF  N 

C 

«NAX=NTEST 

1 00  DO  1 90  N=NZERO,«NAX, 3 

C 

C  SET  F (N) ,F (H-1 ) 

C 

FN1 =1 .OE-28 
FH= .  0 
ALPHA=.0 

IF (N-(N/ 2) *2 >120,11 0,1 20 
110  JT=-1 

GO  TO  130 
120  JT  =  1 
130  H2=H-2 

BO  160  K  -  1  ,M2 
HK-M-K 

BNK  -  2  .*FLQAT (NK)  »FN1 7X-FN 

FN=FM1 

Fill  =BNK 

IF (NK-N-1 >150,140,150 
140  BJ=BHK 
150  JT --JT 
S=1+JT 

160  ALPHA=ALPHA+6HK*S 
BHK =2.*FN1/X-FN 
IF(N)180. 170,180 
170  BJ"BNK 
180  ALPHA=ALPHA+BMK 
BJ=BJ/ALPHA 

IF <ABS(BJ-BPREV>-ABS<D*BJ) >200,200,190 
190  BPREV=BJ 
IER=3 

200  RETURN 
END 


PROGRAM  APBAD  l INPUT , OUTPUT , TAPES- INPUT , TAPE4S0UTPUT ) 
DIMENSION  RCXP ( 7 ) ,RCYP ( 7 ) 

DIMENSION  DXDY (9,1?) 

DIMENSION  EPHIR(5,9) ,EPHII (5,9) 

DIMENSION  EXR(9, 1/),EXI(9, 17),EYR(9, 17)  ,EYI (9,17) 

DIMENSION  RREC(7),RFR(7) 

DIMENSION  SHEF A < 1 00 ) , FREGA( 1 00) 

REAL  HZFRR, HZFRI 

REAL  MCXR (9, 1 7) ,MCXI (9, 17) ,NCYR<9, 17) ,MCYI (9,17) 

REAL  MAG 
REAL  HTHR, HTH1 
COMPLEX  TRAN 
COMPLEX  TRANFR( 9, )7) 

COMPLEX  MX,MY,CEP,FXI,FYI,BARG 
COMPLEX  MCX,MCY,CHZX,CHZY,HZ 
COMPLEX  FLUX,HZF(7> ,HZA (7) 

DIMENSION  WEIT(7) 

REAL  MUO 
COMPLEX  ZHAT 
REAL  LAM 

C0NH0N/BLQCK1/U,F 

COMMON/P ARAMS/’NOPT (7) ,T(7),SIGR(7) ,EMUR(7) ,EPSR(7) ,DLUP(7) , 
+FRCN(7),NG0,LL(7) 

DATA  RFR/. 6096,. 7620,. 9144, 1.0*48,1.2192, 

+1 .8288,2.4384/ 

DATA  RREC/.3048, .4572, .6096, .7620, .91 44, 

♦1 .5240.2.1336/ 

DATA  UEIT/. 25,. 125,. 125,. 125,. 125,. 125,. 125/ 

C  RCXP  AND  RCYP  ARE  THE  POINTS  FOR  INTEGRATION  OVER  THE 
C  RECEIVING  12  IN.  LOOPS 

DATA  RCXP/O.O,. 04205374,. 124 10748,. 06205374, -.04205374, 

+-. 12410748, -.06205374/ 

DATA  RCYP/O. 0,.1 07480231, 0.0,-.t 07480231, -.1074B0231, 0.0, 
+.107480231/ 

PI=3. 1415926536 

TP=2.*PI 

PIQT=PI/2. 

ETA=377. 

MU0=4 .*PI*1 .E-07 

DO  5  1-1,9 

DO  5  J=1 ,17 

TRANFR 1 1 ,  J )=  < 1 . 0,0.0) 

5  CONTINUE 

C  A  IS  THE  AREA  OF  THE  LOOP  ANTENNAS 
A=7.258335668E-2 
C  PROBLEM  LOOP 
KKK=0 

1000  CONTINUE 
KKK=KKK+1 
WRITE  (6,400 )KHK 
CALL  IADATA 
FR=FRCN  < 1 ) 

FRMAX=FRCN(2) 

FRINC=10.**FRCN<3) 

I F  LP =0 
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1234  F=FR*1.E+6 
IFLP=IFLP+1 
FREQA(IFLP)rF 
U=TF*F 

LANr(3.Q£08)/F 
iMAX'NUPT ( I ) 

JrtAX=NQPT ( 2 ) 

DX^BLUP C 1 ) /FLOAT  v IrtAX-1 ) 

DY  =  BLUP(2  )/'FLOAT(  JMAX-f ) 

IMAX21-If1AX/2+1 
IMAX22-I rtAX/2+2 
JrtAX2  J -JHAX/2  +  1 
JNAX22=JHAX/2+2 
C 

C  DIFFERENTIAL  AREA  ASSIGNMENT 

C 

DO  25  1=1 , INAX 
DO  25  J=1 , JNAX 

IF T I  .NE.  IHAX21  .AND.  I  .NE.  IMAX22)  GO  TO  26 
IF< I  .EQ.  INAX21  .OR.  I  .EQ.  IHAX22)  GO  TO  27 

26  IF( J  .EQ.  JNAX21  .OR.  J  .EQ.  JNAX22)  GO  TO  28 
BXDYil, J)=DX*DY 

GO  TO  25 

28  DXDYt I . J  >  =  (DX»DY )/2. 

GO  TO  25 

27  1KJ  .EO.  JNAX2I  .OR.  J  .£0.  JNAX22)  GO  TO  29 
DXDY(  I ,  J )  -  <  I* X * B Y )  /2. 

GO  TO  25 

29  DXDYT l , J>= (DX*DY)/4. 

25  CONTINUE 

L 

C  RFR  IS  THE  ANTENNA  SEPARATION  IN  FREE  SPACE 
C  RREC  IS  THE  DISTANCE  OF  THE  RECEIVING  ANTENNA  FROM  THE  APERTURE 
C 
C 

C  CALCULATE  FREE  SPACE  FLUX  THROUGH  ANTENNA 
C 

DO  20  1=1,1 
DO  10  J=1 ,7 
X=RCXP( J) 

Y=RCYP( J) 

XYS -SORT  <  X=*X+Y  *Y  > 

XYZS=SQRT(X*X+YtY+RFR(I ) **2) 

DEL=SQRT( X  *X+Y*Y ) 

ARG=DEL/RFR( I ) 

TH£TA=ATAN( ARG) 

FFRE=A«COS( THETA) /TP 
R=RFR(I)/COS<THETA) 

BETA=TP/LA« 

BETR=BETA»R 

HTHR=< (-1  .)*<BETA=»*2)+SINUHETA)*A/(4.*PI*R) >* (COS(BETR) 
*-SIN<BETR)/(FETR)-C0S(BE7R)/(BETR**2)) 

HTHI=((-1.)*(BETA*+2 )*SIN(THETA )*A/(4.*PI»R ) )  MSINTBETR  )/(BETR 
+**2)-SIN(BETR)-C0S(BE7R)/(BETR)> 

HZFRR=FFRE*<  COS (BETR )/ <R  **3)+BETA*SIN(BETR »/ (R^*2) )*COS<  THETA) 
HZFRI=FFRE*(BETA*C0S(BETR>/(R**2)-SIN(BETR)/(RM3) >*COS< THETA) 
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10  HZF  ( J )  =CNPLX  <  (HZFRRtRFR  (D/XUS)-  <HTHR*X  YS/XY2S ) , 
♦  <HZFRI+fiFRd;/XYZS>-<HTHI*XYS/XYZS) ) 

FLUX -(0.0, 0-0) 

00  1’J  11-1  r  Z 

15  FlUX=FLUX+A*UE1T  til) *HZF(II) 

HAG-CABStFLUX) 

DB~20 . *ALOG10<MAG) 

20  DB1 =DB 
C 

C  START  TAN  £  FIELD  LOOP  (FIRST  QUAD.  ONLY) 

C 

445  CONTINUE 

DO  50  1=1 , IHAX21 
DO  50  J=1 , JNAX21 
X  =  <  1  - 1  )*DX 
Y-( J-1 ) *DY 
DEL=SDRT<X*X+Y*Y) 

AR6=Dt'L/(  .305 ) 

THETA=ATAN(ARG) 

CQTHE-COS(THETA) 

SITHE=SIN( THETA) 

R=(.305)/CQTHE 
BE  IR»(TP:*R)/IAN 
COSBR=COS(BETR ) 

SINBR=SIN(BETR) 

FACT- (ETA*PI*A*SI THE )/(LAH*LAH*R) 

EPHIRd,J)  =  FACT*<COSBR-(SINBRXBETR)) 

EPHIId,  J)=(-1.)+FACT*(SINBR+(C03BR/BETR)) 
IF(EPHIR<I, J)  .LT.  1.0E-05)  EPHIR( I , J)=0.0 
IF(I.EQ.1)G0  TO  40 
PHI=ATAN(Y/X) 

40  IFU.EQ.1  )PHI=PIQT 

IF (NOPT (4)  .EQ.  0)  60  TO  456 

U=(TP/LAM)*SITHE*COS(PHI) 

U=(TP/LAN)*SITHE*SIN(PHI) 

CALL  SHIELB(U,V,R,TRAN) 

TRANFRd,  J)=TRAN 

510  F0RNAT(/,2X, * TRAN =*,2E1 4.7) 

456  CONTINUE 

EXR(I  ,J)=EPHIR(I  ,J)+SIN(PHI) 

EXI(I,J) -EPHIId  jJXSINtPHI) 

EYRd , J)=EPHIR(I,  J)*COS(PHI) 

EYI (I , J) =EPHII(I , J)*COS(PHI ) 

IF (I  .EQ.  1)  EYR(I,J)=EYId, J)=0.0 
50  CONTINUE 
C 

701  FORNAT  (/,2Xf:*FLAG1*) 

C  END  TAN  E  FIELD  LOOP 
C 

C  FILL  UP  OTHER  QUADRANTS 
C 

DO  1  J=1 ,  JNAX21 
DO  1  1-2 , 1 MAX  2 1 
II=INAX-(I-2) 

TRANFR(II,J)=TRANFR(I, J) 

EXR(II , J)rEXR(I, J) 

EXI(II,J)=EXId,  J) 

EYR( II ,  J )-EYR(I , J) 
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n  n  n 


1  EYKI  I ,  J>=ETf I  (I .  J» 

DO  2  1=1 , 1  MAX  2 1 
DO  2  J=2,J«AX21 
JJ=JNAX-vJ-2) 

TRANFRU,JJ)=TRANFR(I,J> 

E  XR ( I ,  JJ)=EXR\ I, J ) 

EXI ( I , J J ) =EX I ( I ,  J  i 
EYR(I,JJ)— EYR(I rJ) 

2  £ r I < I , JJ )  --EYKI, J) 

DO  3  1=2, INAX21 
DO  3  J=2, JNAX21 
1 1  =  I MAX2 1  -*■  I 
JJ= JNAX21 +J 

TRANFR1 II , JJ)=TRANFR( I , J) 

E X ft  1 1 I ,  JJ)=EXR(I , J ) 

EXI < II, JJ)=EXI ( I,  J) 

EYR < II , JJ)=EYR(  I , J ) 

3  EYI  i  II , JJ)=EY1 ( I , J) 

CORRECT  SIGNS  OF  X  S  Y  COrtPQNENTS 

DO  60  1  =  1 ,  IHAX 
DO  60  J=1 , JNAX 

IF ( I  .LT.  INAX22  .AND.  J  .LT.  JNAX22)  GO  TO  61 

IF ( I  .LT.  INAX22  .AND.  J  .6E.  JNAX22)  60  TO  60 

IF ( I  .GE.  INAX22  .AND.  J  .IT.  JHAX22)  GO  TO  62 

IF  \  I  .GE.  INAX22  .AND.  J  .GE.  JNAX22 )  GO  TO  63 

61  EXR<I,J)  =  (-1.)*EXRU,J) 

EXI(I,J)=(-1. )*EXI( I, J) 

GO  TO  60 

62  EXR(I,J)=<-1 . )»EXR(I,J) 

EXI(I,J)=i-t.)*EXI(I,J) 

EYR(I,J)=(-1.)»EYR(I,J) 

EYI(I,J)  =  1-1.)*EYK1,J> 

GO  TO  60 

63  EYR(I,J)=(-I.)*EYR(I,J) 

EYi(I,J)=(“1.)*EYI(I,J) 

60  CONTINUE 
C 

C  CALCULATE  NAGNETIC  CURRENT 

f. 

DO  70  1=1 , INAX 
DO  70  J=1 , JNAX 
HCXR(I,J)=2.*EYR(I,J) 

MCXI (I ,  J)=2.*EYJ  <1 ,  J) 

HCYR( I ,J)  =  <-2.  )+EXR< I ,J) 
MCYIU,J)=(-2.)*EXI<I,J) 

70  CONTINUE 
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C  LALCU1.AIL  l  COMPONENT  OF  H  FIELD  At  HUINI  XPI',fPP  IN  RECEIVER  H-ANt 

<; 

DO  80  K=1 , 1 
HO  90  IR=1 ,7 
XPP=RCXP( IR) 

YPP=RCYP( IR) 

HZA<IR>=CHPLX<0. 0,0.0) 

DO  30  1=1 ,  IHAX 
DO  30  J=1 ,  JHAX 
C 

C  APERTURE  COORDINATES 
C 

IF C I  .LT.  IHAX22  .AND.  J  .LT.  JHAX22)  GO  TO  31 

IF C I  .LT.  IMAX22  .AND.  J  .GE.  JNAX22)  GO  TO  32 

IF ( I  .GE.  IHAX22  .AND.  J  .GE.  JHAX22)  GO  TO  33 

IF< I  .GE.  IHAX22  .AND.  J  .LT.  JNAX22)  GO  TO  34 

31  XP=<I-1)*DX 
YP=( J*1 )*DY 
GO  TO  35 

32  XP=(I-1)-»DX 

YP=  <  —  1 . )*( JHAX-J+1 )*DY 
GO  TO  35 

33  XP=(-1  .)••»(  INAX-I  +  1  )*DX 
YP=M.)*(JNAX-J+1)*DY 
GO  TO  35 

34  XP=(-1 . )*( INAX-I+1 )*DX 
YP= ( J-1 )*DY 

35  RHRP=SQRT  ( <  XP-XPP )  *  *2  +  ( YP-YPP )  *=*2+RREC ( K )  r*2 ) 

DENQH=RHRP**5 

TK=(TP/LAH)»RHRP*(-1 .0) 

BARG=C«PLX(0.0,TK) 

CEP-CEXP ( BARG ) 

HCX=CHPLX<NCXR(I  ,J) ,NCXI U , J) )*TRANFR< I , J) 

HCY=CHPLX(HCYR(I, J) , MCYI< I, J) )*TRANFR< I , J > 

C 

C  CHZX  IS  THE  CONTRIBUTION  TO  HZ  FROM  KCX  tX  CONP.  OF  NAG.  CURRENT) 

C  CHZY  IS  THE  CONTRIBUTION  TO  HZ  FROH  HCY  (Y  CONP.  OF  HAG.  CURRENT) 

C 

CHZY=3.*(YPP-YP)*RREC(K)*NCY*CEP*DXDY(I,J)/(DEN0N*4.*FI) 

CHZX=3.*(XPP-XP)*RREC(K)*HCX*CEP*DXDYU,J)/(DENQN*4.*PI) 

ZHAT*CHPLX<0.0,(-1.)/(U*HU0)) 

CHZX=CHZX*ZHAT 

CHZY=CHZY*ZHAT 

30  HZA( IR)=HZA( IR)+CHZX+CHZY 
90  CONTINUE 

FLUX=<0. 0,0.0) 

DO  95  JI=1,7 

95  FLUX=FLUX+A*UEIT<JI)*H2AUI) 

HA6=CABS(FLUX) 

DB=20.*AL0G10<HAG) 
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80  062-flB 

SHEF-DB1  -DB2 
SHEFA! IFLP)=SHEF 
FR=FR*FRINC 

IF (Fft  .LE,  FRNAX )  GO  TO  1 23-4 
URITEI6 ,104) 

URITE(6,105)(FREQA(I) ,  SHEFA < I ) ,  1-1  ,IFLP) 

444  CONTINUE 

IF  (NGO  .EO.  1)  GO  TO  1000 
100  FORMAT </,  15X ,*X  COMPONENT*) 

200  FORMAT (/, 15X ,*Y  COMPONENT*) 

209  FORMAT </,2X,*P0INT*,F5.2r1X,F5.2) 

210  FORMAT ( / , 2X , *DISTANCE=*,F7 .4 ,2X, *HZ~* ,2E14.6 ) 

300  FORMAT (/,2X.*PHI  COMPONENT*) 

301  FORMAT!/, 2X, *1  * , 3X, *J* ,3X  ,*DEL* ,5X, *THE7A*> 

302  F0RHAT(/,2X, I2,3X, I2,3X,F?.4,5X,F7.4) 

303  FORMAT l/,2X, +PHI  CQMP.  INSIDE  LOOP*) 

304  F0RMAT(/,2X,I2,3X,I2,3X,EI4.7,2X,E14.7) 

305  FGRHAT(/',2X,I2,3X,I2,3X,*PHI=*,F8.4,*PI*) 

306  FORMAT!/, 3X,*MAG.  CURRENT*) 

30?  FORMAT (/, 3X, *1 *, 3X, * J*,3X,*MCXK*,14X, *NCXI*, l 4X,*MCYR*,14X, *MCYI») 

308  FORMAT !/,3X,I2,3X, 12,4! El  4.?)) 

309  F0RMAT!/,3X,I2,3X,I2,4!EM.?)> 

310  FORMAT </,3X,*I*,3X,*J*,3X,*EXR*,14X,*EXI*,14X,*EYR*,14X, *EYI*) 

311  FORMAT!/, 2X,*FREE  SPACE*) 

312  FORMAT!/, 2X , *APERTUR£* ! 

313  FORMAT !  /  ,2X,;*D  I  STANCE—*,  F7.4,2X,*HZFR=*  ,2E  M.6) 

314  F0RMAT(/,2X,*DISTANCE=*,F7.4,2X,*FLUX=*,2E14.6,2X, »DB=*,F10,4) 

400  FORMATdHI ,/////, 30X , fPROBLEM  *,I2,*  INPUT  DATA  *,//> 

500  FORMAT!/, 2X,*HELPt) 

501  FORMAT!/, 2X,*HELP1*) 

104  FORNAT ( 1  HI, 30X,*FREQ<HZ)*, 1 0X, *SHEF«DB )♦,//) 

105  F0RMAT(30X,E10.3,5X,F10.3) 

999  CONTINUE 

STOP 

END 

SUBROUTINE  SHIELD(U,V,R,TRAN) 

COMMON/BLOCK 1 /U ,F 

COMMON/PARAi1S/NOPT!7) ,T(?) ,SIGR(?) ,EMUR( ? ) ,EPSR(?) , 

+DLUF (7) ,FRCN<?) ,NGO,LL<?) 

COMPLEX  CEP.CEN 
COMPLEX  TOP.BOI 
COMPLEX  TRAN, PROD 

COMPLEX  GS(7),TAU(7),PSY(7),HT(7),SZR(7),P0LE(7),CSPL(7) 

COMPLEX  2(7) 

REAL  LAM 

PI  =  3. 1415926536 

EMUOd  .26E-6 

EPS0=8.854E-12 

SIGC=5.80E+7 

LAM -!3.0E08)/F 

RO-R 
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_ (1C  IS  MIL  IUIAI  NO.  OF  LAYERS 

HC^NQP  T  v3 ) 

MCH--HC-1 

- CALCULATE  LAYER  IMPEDANCES 

DO  5  M=1f MC 

EMU=EMUR(M)*EMUO 

EPS=EPSR<M) +EPSO 

S1G=SIGR(H) *SIGC 

IF CSIGft(M)  .NE.  O.O)  GO  TO  6 

B£TA=(2.*PI)/LAN 

IF (MC  .EQ.  5  .AND.  M  .EO.  3)R=R0+T (2)+T(3) 

IF1MC  .EO.  5  .AND.  M  .EQ.  5 ) R-RO+T { 2 ) +T (3)+T(4) 

IF (MC  .EQ.  7  .AND.  M  .EQ.  3)R=R0+T(2)+T(3) 

IF ( MC  .EQ.  7  .AND.  H  .EQ.  5)R=R0+T(2>+T(3)+T(A)+T(5) 

IF(MC  .EQ.  7  .AND.  H  .EQ.  7)R=R0+T(2)+T(3)+T(4)+T(5)+T(6) 

TOP  ~CMPL  X ( (-1.0)*<BETA**2)*1R**2i ,BETA»R) 

B0T  =  CMPLXU.-(BETA**2)*(R**2) ,  BET  A*R ) 

2 ( M) = i 377. ) ♦ ( TOP/BOT ) 

100  F ORMAT (/ ,2X, I2,5X, 2E1 4.6) 

GO  TO  5 

6  SKDP=( 1 .O./iSQRT  <PI*F*EHU*SIG) ) 

Z  (M)~CMPLX<  I .  /  (SIG»SKDP) ,  1 ./  (SIG*SKDP) ) 

5  CONTINUE 
DO  10  M=1,MC 
EMU=EMUR(N)*EHU0 
EPS=EPSR ( H )*EPSO 
SIG=SIGR(N)*SIGC 
GNR=-U*W*EI1U*EP3 
GH1-U*EHU*SIG 

10  GS<N)=CNPLX(GMR,GMI) 

AL=SQRT(U*U+V*V.' 

HT<MC)=0. 

S2R<HC)-I. 

DO  11  (1=1,  MC 

11  TALK  M )=CSQRT (GS(M)+U»U+V»V) 

DO  12  M--2,NCM 
PaY<i1)-TAlj(M/ *  I  i (1  > 
cep-ceaF'  (F"jY  in>) 

CEM-1 ./CEP 
CSPLu1)='CEP+CEM)/2. 

HT(M)=(CEP-CE«)/(CEP+CEM) 

12  CONTINUE 

DO  13  HR=t,MCM 

M--rtCi1-(1R+1 

HP-H+1 

13  SZR<M)=(2{MF)/Z(N))*<SZR(MP)+HT(MP))/( 1 ,+SZR (MP)*HT(HP) ) 

POLE  < 1 )  =  1 .+SZR1 1 ) 

PR0D=2.t(Z<NC)/Z(  1 )  F/POI.E  ( 1 ) 

DO  14  M=2,NCM 

POL  E  <  N ) 3 1 .4SZR(M(*HI<M) 

14  PROO=PROD/(POLE(M)*CSPt <H )  > 

TRAN=PROD 

RETURN 

END 

SUBROUTINE  IADATA 

COMMON/PARAMS^NOPT (7) ,T(7) ,SIGR(7) ,ENUR(7 ) , EPSR (7) ,DLUP(7) ,FRCN(7) 
1 ,NG0,LL<7) 
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NJO-y 

WRITE  (6,600) 

READ  (5,500)FRCN( 1 ) 

URITE (6,607) 

R£AD(5,5Q0>FRCN(2) 

FRCN(3)=. 125 
URITE  (6,60t ) 

READ  (5,501 )DLUP(1 ) 

URITE  (6,651) 

READ  (5,501 )DLUP(2) 

URITE  (6,602) 

READ  (5.502)N0P7(1) 

URITE  (6,652) 

READ  (5,502)N0PT (2) 

URITE  (6,60J ) 

READ  (5,503)NQPT(4) 

IF(NOPT< 4)  ,EQ.  0)  BO  TO  30 
URITE  ( 6,604) 

READ  <5,504)N0PT (3) 

MC=NOPT (3) 

MCM=HC-1 
DO  10  1=2, NCN 
URITE  (6,605)1 
READ  (5,505)7(1) 

READ  (5,505)SIGR(I) 

10  CONTINUE 
30  CONTINUE 

T  (1  )=T(NC)r-.305 
S1GR(1 )=5IGR(MC)=0.0 
URITE  (6,606) 

READ  (5,506 )NGO 
DO  20  1  =  1, HC 
ENUR( I >  =  1 .0 
20  EPSR(  D  -1 .0 

600  FORMAT i/,2X,*EN7ER  START  FREQ  (MHZ)*) 

601  F ORH AT (/, 2 X,* APERTURE  HEIGHT*) 

651  FORMAT (/,2X,*APER7URE  UIDTH+) 

602  FORMAT </,2X,*X  POINTS*) 

652  F0RHAT(/,2X,*Y  POINTS*) 

603  FORMAT!/ ,2X,*SHIELDING  NATERIAL  IN  APERTURE?  Y=1,N=0*) 

601  FORMAT!/, 2X,*H0U  MANY  LAYERS?  (TYRAS  M+1)*) 

605  FORMAT!/, 2X,*ENTER  THICKNESS  AND  SIBR  OF  LAYER  NO.*, IX, 12) 

606  FORMAT'/, 2X,*M0RE  DATA  LATER?  Y  =  1 ,  N=0*) 

607  FORMAT!/, 2X,*HAX  FREQ  <NHZ>*> 

500  FORMAT(E)O.I) 

501  FORMAT (2E1 0.  i ) 

502  FORMAT (12) 

503  F  ORM AT (12) 

504  FORMAT! 12) 

505  FORMAT (F14.8) 

506  FORMAT (12) 

RETURN 

END 
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